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Abstract 



We present a charge-dependent nucleon-nucleon (NN) potential that fits the 
world proton-proton data below 350 MeV available in the year of 2000 with 
a x 2 P er datum of 1.01 for 2932 data and the corresponding neutron-proton 
data with \ 2 /datum = 1.02 for 3058 data. This reproduction of the NN 
data is more accurate than by any phase-shift analysis and any other NN 
potential. The charge-dependence of the present potential (that has been 
dubbed 'CD-Bonn') is based upon the predictions by the Bonn Full Model 
for charge-symmetry and charge-independence breaking in all partial waves 
with J < 4. The potential is represented in terms of the covariant Feynman 
amplitudes for one-boson exchange which are nonlocal. Therefore, the off-shell 
behavior of the CD-Bonn potential differs in a characteristic and well-founded 
way from commonly used local potentials and leads to larger binding energies 
in nuclear few- and many-body systems, where underbinding is a persistent 
problem. 

PACS numbers: 13.75.Cs, 21.30.Cb, 25.40.Cm, 25.40.Dn, 24.80.+y 
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I. INTRODUCTION 



In the 1970's and 80's, a comprehensive fieldtheoretic meson-exchange model for the 
nucleon-nucleon (NN) interaction was developed at the University of Bonn. The final 
version, published in 1987, has become known as the Bonn Full Model 0. For a pedagogical 
review see Ref. 0. 

In the language of fieldtheoretic perturbation theory, the lowest order contributions to 
the NN interaction generated by mesons are the one-boson exchange diagrams. Further- 
more, there are many irreducible multi-meson exchanges. The diagrams of 2tt exchange are 
most prominent since they provide the intermediate-range attraction of the nuclear force. 
However, once explicit diagrams of 2tt exchange (with intermediate A isobars) are used in 
a model, then it is vital to also include the corresponding diagrams of up exchange. There 
are characteristic (partial) cancellations between the two groups of diagrams, which are 
crucial for a quantitative reproduction of the NN data. Moreover, the Bonn model con- 
tains additional classes of irreducible 3tc and 4tc exchanges which are important conceptually 
rather than quantitatively, since they appear to indicate convergence of the diagrammatic 
expansion chosen by the Bonn group [|TJ. 

The development of the Bonn Full Model was necessary to test reliably the meson- 
exchange concept for nuclear forces and to assess systematically the range of its validity. 
Thus, the model represents a benchmark for any alternative attempt (based, e. g., on quark 
models, chiral perturbation theory, or other ideas) to explain the nuclear force. 

Due to its comprehensive character, the Bonn model provides a sound basis for addressing 
many important issues. One of them is the charge dependence of nuclear forces. The charge- 
symmetry-breaking (CSB) of the NN interaction due to nucleon mass splitting has been 
investigated in Ref. . It turns out that considerable CSB is generated by the 27r-exchange 
contribution to the NN interaction and the irp diagrams such that the CSB difference in 
the singlet scattering lengths can be fully explained from nucleon mass splitting. Also, 
noticeable CSB effects occur in P and D waves. Empirical evidence for CSB is seen in the 
the Nolen-Schiffer (NS) anomaly [|J regarding the energies of neighboring mirror nuclei. A 
recent study || has shown that the CSB in partial waves with L > as derived from the 
Bonn model is crucial for a quantitative explanation of the NS anomaly. 

The charge-independence breaking (CIB) of the NN interaction has also been inves- 
tigated |J. Pion mass splitting is the major cause, and it is well-known that the one- 
pion-exchange (OPE) explains about 50% of the CIB difference in the singlet scattering 
lengths. However, the 27r-exchange model and the diagrams of three and four irreducible 
pion exchanges contribute additional CIB which can amount up to 50% of the OPE CIB- 
contribution, in S, P, and D waves. This effect is not negligible. 

Other important issues related to the nuclear force are relativistic effects, medium ef- 
fects, and many-body forces to be expected in the nuclear many-body problem. The medium 
effects on the nuclear force when inserted into nuclear matter have been calculated thor- 
oughly. A large repulsive contribution to these medium effects comes from intermediate A 
isobar states which also give rise to energy-dependence. On the other hand, isobars create 
many-body forces that are attractive. Thus, large cancellations between these two classes 
of many-body forces/effects occur and it has been shown that the net contribution is very 
small 0. Relativistic effects, however, may play an important role in the nuclear many-body 
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problem [0. 

Multi-meson exchange diagrams are very involved. Moreover, contributions of this kind 
are, in general, energy dependent. This would make the NN potenital — defined as the sum 
of irreducible diagrams — energy dependent. A NN potential that depends on energy creates 
conceptual and practical problems when applied in nuclear many-body systems. For a large 
class of nuclear structure problems, these complications are without merit. 

For these reasons, already early in the history of the meson theory of nuclear forces, the 
so-called one-boson-exchange (OBE) model was designed which — by definition — includes 
only single-meson exchanges (which can be represented in an energy- independent way). 
Usually, the model includes all mesons with masses below the nucleon mass, i. e., tt, rj, 
p(770), and u(782) |§. In addition, the OBE model typically introduces a scalar, isoscalar 
boson — commonly denoted by o (or e). Based upon what we discussed above concerning 
multi-meson exchange contributions, it is clear now that this a must approximate more 
than just the 2tt exchange. In particular, it has to simulate 2tt + irp exchanges which are 
clearly not of purely scalar, isoscalar nature. Consequently, the a approximation is poor (as 
demonstrated in Fig. 11 of Ref. Jl|). One way to make up for this deficiency is to readjust 
the parameters of the a boson in each partial wave. Moreover, the 2tt + irp exchanges 
create — in terms of ranges — a very broad contribution that cannot be reproduced well by a 
single boson-mass; two masses will do better. The fact that we are dealing here with a very 
broad mass distribution is supported by an entry in the Particle Data Tables || which lists 
a cr (or fo) with a mass between 400 and 1200 MeV. 

Based upon the philosophy just outlined, we have constructed a NN potential that is 
energy-independent and defined in the framework of the usual (nonrelativistic) Lippmann- 
Schwinger equation. Thus, it can be applied in the same way as any other conventional 
NN potential. The crucial point, however, is that it reproduces important predictions by 
the Bonn Full Model, while avoiding the problems that the Bonn Full Model creates in 
applications. The charge-dependence (CD) predicted by the Bonn Full Model is reproduced 
accurately by the new potential, which is why we have dubbed it the CD-Bonn potential. 
The off-shell behavior of CD-Bonn is based upon the relativistic Feynman amplitudes for 
meson-exchange. Therefore, the CD-Bonn potential differs off-shell from conventional NN 
potentials — a fact that has attractive consequences in nuclear structure applications. 

An earlier version of the CD-Bonn potential — which, however, did not contain all the 
charge- dependence — was published in Ref. |f§ where the off-shell aspects are discussed in 
great detail. 

In Sec. II, we present the potential model. Charge dependence is discussed in Sec. III. 
The results for NN scattering and the deuteron are presented in Sec. IV and V, respectively. 
Conclusions are given in Sec. VI. The paper has four appendices which contain mathematical 
and other details. 

II. THE MODEL 

As discussed in the Introduction, the CD-Bonn potential is based upon meson exchange. 
We include all mesons with masses below the nucleon mass, i. e., tt, rj, p(770), and c<j(782). 
Besides this, we introduce two scalar-isoscalar a bosons. 
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For the rj (with a mass of 547.3 MeV), we assume a vanishing coupling to the nu- 
cleon, which implies that — de facto — we drop the 77. This assumption is supported by 
semi-empirical evidence from various sources. Analyzing NN scattering data in terms of 
forward dispersion relations, Grein and Kroll |l(J determined the t]NN coupling constant 



to be consistent with zero. Tiator and coworkers [I1J extracted the 77 coupling from rj pho- 



toproduction data and found <?~/47T = 0.4. Such a small coupling constant generates a 



negligible contribution in the NN system. In the development of the Bonn Full Model for 
the NN interaction [|lj, it was noticed that a good fit of the NN data favors a vanishing rj 
contribution. 

In Table [I], we list the hadrons involved in our model together with their masses and 
coupling parameters. For the nNN coupling constant, we choose the 'small' value S^/47T = 



13.6 — consistent with recent determinations by the Nijmegen JT^|13| and VPI group fI~4}-|T6"fl . 
It is appropriate to mention that the precise value of the ttNN coupling constant is an 
unsettled issue at this time, and we refer the interested reader to Refs. [0,0] f° r a critical 
discussion and review of the topic. For the vector mesons p and u, for which precise empirical 
determinations of the coupling constants are difficult (if not impossible), we use the values 
from the Bonn Full Model 0J. 

We start from the following Lagrangians that describe the coupling of the mesons of 
interest to nucleons: 

£ir°NN = -g^i^r^if^ (2.1) 

C W ± NN = -V2 g7r ±^ 5 T ± ^ {7T±) (2.2) 

C.nn = -gM^ (2.3) 

CuNN = -g^Y^ ] (2.4) 



C pNN = -g p W>T$ ■ tpf - J^f^^ ■ (9M P) ~ b&f) (2-5) 

where ip denotes nucleon fields, (p meson fields, and t^± are standard definitions of Pauli 
matrices and combinations thereof for isospin | M p is the proton mass which is used 



as scaling mass in the pNN Lagrangian to make f p dimensionless. To avoid the creation of 
unmotivated charge dependence, the scaling mass M p is used in the pNN vertex no matter 
what nucleons are involved. 

In the c. m. system of the two interacting nucleons, the OBE Feynman amplitude gen- 
erated by meson a is, 

- iv a u, q ) = ^(qo^^i/^^c-qo^^c-q) (2 6) 

(q' - q) 2 - m l a 

where (i = 1,2) are vertices derived from the above Lagrangians, Ui Dirac spinors 
representing the interacting nucleons, and q and q' their relative four-momenta in the initial 
and final states, respectively; P a divided by the denominator is the appropriate meson 
propagator. 

The one-boson-exchange potential is defined by (i times) the sum over the OBE Feynman 
amplitudes of the mesons included in the model (Fig. [I]); i. e., 

V(q!, q) = MM E V a {q[, q)^(q', q; A Q ) . (2.7) 

O=7T jTT^ ,p,W,(Tl ,(T2 
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As customary, we include form factors, J- a (q!, q; A a ), applied to the meson-nucleon vertices, 



and a square-root factor M/\/E'E (with E = y/M 2 + q 2 , E' = \J M 2 + q' 2 , and M the 
nucleon mass). The form factors [see Appendix B, Eq. (|B9|), for details] regularize the 
amplitudes for large momenta (short distances) and account for the extended structure 
of nucleons in a phenomenological way. The square root factors make it possible to cast 
the unitarizing, relativistic, three-dimensional Blankenbecler-Sugar (BbS) equation |20| for 
the scattering amplitude [a reduced version of the four-dimensional Bethe-Salpeter (BS) 



equation [21|] into the following form (see Appendix A for a proper derivation): 

T(q', q) = V(q[, q) + J d 3 k V(q', k) 2 _^ 2+ . £ T(k, q) (2.8) 

Notice that this is the familiar (non-relativistic) Lippmann-Schwinger equation. Thus, 
Eq. (|2.7| ) defines a relativistic potential which can be consistently applied in conventional, 
non-relativistic nuclear structure, in the usual way. 



The Feynman amplitudes, Eq. (2.6), are in general nonlocal expressions; i. e., Fourier 
transform into configuration space will yield functions of r and r', the relative distances 
between the two in- and out-going nucleons, respectively. The square root factors in Eq. ( |2.7|) 
create additional non-locality. 

While for heavy vector-meson exchange (corresponding to short distances) non-locality 
appears quite plausible, we have to stress here that even the one-pion-exchange (OPE) 
Feynman amplitude is non-local. This fact is often overlooked. It is important because 
the pion creates the dominant part of the tensor force which plays a crucial role in nuclear 
structure. 

Applying the nNN Lagrangian Eq. Q2.1|) to the amplitude Eq. (|2.6| ) yields the one- 
pion-exchange (OPE) potential (suppressing charge- dependence and isospin factors for the 
moment) 

v , / )= 9l (E' + M)(E + M) ( gW _ <Ti-q \ 
Aq ' qj AM 2 (q' -q) 2 + ml \E' + M E + MJ 

\E' + M E + MJ K ' 

If we would now apply the approximation, E' ~ E ~ M (static approximation), then this 
simplifies to 

v (ioc ](k) _ 9l (<n-k)(cr a -k) 

W_ AM 2 k 2 + m 2 n (2 ' iUj 

with k = q' — q. Fourier transform of this latter expression yields, 




.,, • (2.11) 

m n r [m^r) 2 I r 
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This is the local OPE potential that is used by most practitioners. However, the important 
point to notice here is that this local OPE is not the full, original OPE Feynman amplitude; 
it is an approximation. 

The obvious question to raise at this point is: How much does the local approximation 
change the original result or, in other words, how drastic is the local approximation? 

For this purpose, we show in Fig. ^| the half off-shell 3 S\- 3 Di potential that can be pro- 
duced only by tensor forces. The on-shell momentum q' is held fixed at 265 MeV (equivalent 
to 150 MeV lab. energy), while the off-shell momentum q runs from zero to 2000 MeV. The 
on-shell point (q = 265 MeV) is marked by a solid dot. The solid curve is the relativistic 
OBE amplitude of it + p exchange. Now, when the relativistic OPE amplitude, Eq. ( |2.9| ), is 
replaced by the static/local approximation, Eq. Q2.10D , the dashed curve is obtained. When 
this approximation is also used for the one-p exchange, the dotted curve results. It is clearly 
seen that the static/local approximation does change the potential drastically off-shell: it 
makes the tensor force substantially stronger off-shell. 

In summary, one characteristic point of the CD-Bonn potential is that it uses the Feyn- 
man amplitudes of meson exchange in its original form; local approximations are not applied. 
This has impact on the off-shell behavior of the potential, particularly, the off-shell tensor 
potential. It is well known that the off-shell behavior of an iViV potential is an impor- 
tant factor in microscopic nuclear structure calculations. Therefore, the predictions by the 
CD-Bonn potential for nuclear structure problems differ in a charcteristic way from the 
ones obtained with local NN potentials. For more discussion of this issue, see Sec. VI and 
Refs. |M. 



III. CHARGE DEPENDENCE 



By definition, charge independence is invariance under any rotation in isospin space. 
A violation of this symmetry is referred to as charge dependence or charge independence 
breaking (CIB). Charge symmetry is invariance under a rotation by 180° about the ?/-axis in 
isospin space if the positive z-direction is associated with the positive charge. The violation 
of this symmetry is known as charge symmetry breaking (CSB). Obviously, CSB is a special 
case of charge dependence. 

CIB of the strong NN interaction means that, in the isospin T = 1 state, the proton- 
proton (T z = +1), neutron-proton (T z = 0), or neutron-neutron (T z = —1) interactions are 
(slightly) different, after electromagnetic effects have been removed. CSB of the NN interac- 
tion refers to a difference between proton-proton (pp) and neutron-neutron (nn) interactions, 
only. For recent reviews on these matters, see Refs. fl2"3"|j2"4| . 

CIB is seen most clearly in the 1 Sq NN scattering lengths. The latest empirical values 
for the singlet scattering length a and effective range r are: 



i N 
pp 

,N 



"np 



-17.3 ±0.4 fm [24], 

-18.9 ±0.4 fm mm, 

-23.740 ± 0.020 fm IE71 



r N 
pp 



np 



2.85 ± 0.04 fm |J 



2.75 ±0.11 fm 0; 
2.77 ±0.05 fm |27|. 



(3.1) 
(3.2) 
(3.3) 



The values given for pp and nn scattering refer to the nuclear part of the interaction as 
indicated by the superscript A; i. e., electromagnetic effects have been removed from the 
experimental values. 
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The above values imply that charge-symmetry is broken by the following amounts, 




(3.4) 
(3.5) 



AdciB = dpp — d; 

Ar C iB =rl- r. 



np 



np 



6.44 ± 0.40 fm, 
0.08 ± 0.06 fm. 



(3.6) 
(3.7) 



In summary, the NN singlet scattering lengths show a small amount of CSB and a clear 
signature of CIB. 

The current understanding is that — on a fundamental level — the charge dependence of 
nuclear forces is due to a difference between the up and down quark masses and electromag- 
netic interactions among the quarks. As a consequence of this — on the hadronic level — major 
causes of CIB are mass differences between hadrons of the same isospin multiplet, meson 
mixing, and irreducible meson-photon exchanges. 



The difference between the masses of neutron and proton represents the most basic cause 
for CSB of the nuclear force. Therefore, it is important to have a very thorough accounting 
of this effect. 

The most trivial consequence of nucleon mass splitting is a difference in the kinetic 
energies: for the heavier neutrons, the kinetic energy is smaller than for protons. This 
raises the magnitude of the nn scattering length by 0.26 fm as compare to pp. The nucleon 
mass difference also affects the OBE diagrams, Fig. [I], but only by a negligible amount. 
In summary, the two most obvious and trivial CSB effects explain only about 15% of the 
empirical Aclcsb (cf. Table [Tip. Usual models for the nuclear force include only the two CSB 
effects just dicussed and, therefore, do not reproduce the empirical CSB. 

However, in Ref. || it was found that the irreducible diagrams of two-boson exchange 
(TBE) create a much larger CSB effect than the OBE diagrams and, in fact, fully explain 
the empirical CSB splitting of the singlet scattering length. The major part of the CSB 
effect comes from diagrams of 2n exchange where those with iVA intermediate states make 
the largest contribution. The CSB effect from irreducible diagrams that exchange a 7r and 
p meson were also included in the study. The np diagrams give rise to non-negligible CSB 
contributions that are typically smaller and of opposite sign as compared to the 2-7T effects. 
The net effect explains AacsB quantitatively. 

The above-mentioned investigation || was based upon the Bonn Full Model [Q]. This 
model uses the irNN coupling constant g^/47r = 14.4 which is not current. For that reason 
we have revised the Bonn Full Model using g^/^ti = 13.6 and then repeated the CSB 
calculations of Ref. ||. The total AacsB predicted by the revised model is 1.508 fm (about 
5% less than what was obtained in Ref. [Q with the original model), implying a TBE effect 
of 1.275 fm. 

The only reliable empirical information about CSB of the NN interaction is the scattering 
length difference in the 1 S state, Eq. ( p.4|) . As discussed, the TBE model of Refs [l],[J can 



A. Charge symmetry breaking 



S 



explain this entirely from nucleon mass splitting. For this reason, we have confidence in the 
CSB predictions by this model. Therefore, we will use its predictions also for energies and 
states where no empirical information is available; namely, higher energies in the 1 Sq state 
and partial waves other than 1 So. 

Thus, using the revised Bonn Full Model, we have calculated the difference nn phase shift 
minus pp phase shift without electromagnetic interactions, 5 nn — 5 PP , that is caused by CSB 
of the strong nuclear force due to nucleon mass splitting. The total effect obtained is listed 
in the last colum ('Total') of Table |H| for energies up to 300 MeV and partial wave states in 
which these effects are non-negligible. In that table, we also list the very small effects from 
the OBE diagrams (Fig. [I]) and the kinematical effects (column 'Kinematics') ||30|| . CSB 
phase shift differences are plotted in Fig. |||. It is clearly seen that in most states the TBE 
effect is the largest and, therefore, certainly not negligible as compared to the other CSB 
effects shown. 

Because of the outstanding importance of the CSB effect from TBE, we include it in our 
model. By doing so, we go beyond what is usually done in charge- dependent iViV potentials. 
In most recent models, only the kinematical effects and the effect of nucleon mass splitting 
on the OBE diagrams are included. However, as discussed, this does not explain the CSB 
scattering length difference. Thus, some models leave CSB simply unexplained [J31|, while 



other models add a purely phenomenological term to the potential that fits Aclcsb |32| . 

Before finishing this subsection, a word is in place concerning other mechanisms that 
cause CSB of the nuclear force. Traditionally, it was believed that p° — u mixing explains 
essentially all CSB in the nuclear force |[24|| . However, recently some doubt has been cast on 
this paradigm. Some researchers p3|-|36|1 found that p° — u exchange may have a substantial 
q 2 dependence such as to cause this contribution to nearly vanish in NN. Our finding that 
the empirically known CSB in the nuclear force can be explained solely from nucleon mass 
splitting (leaving essentially no room for additional CSB contributions from p° — u mixing 
or other sources) fits well into this scenario. On the other hand, Miller [23| and Coon and 
coworkers [37| have advanced counter-arguments that would restore the traditional role of 
p-uj exchange. The issue is unresolved. Good summaries of the controversial points of view 
can be found in Refs. p3| , |38|j39|] . We do not include p — to mixing in our model. 

Finally, for reasons of completeness, we mention that irreducible diagrams of tc and 7 
exchange between two nucleons create a charge-dependent nuclear force. Recently, these 
contributions have been calculated to leading order in chiral perturbation theory fjOfl . It 
turns out that to this order the 7r7 force is charge-symmetric (but does break charge inde- 
pendence). 



B. Charge independence breaking 

The major cause of CIB in the NN interaction is pion mass splitting. Based upon the 
Bonn Full Model for the NN interaction, the CIB due to pion mass splitting has been 
calculated carefully and systematically in Ref. ||. 

The largest CIB effect comes from the OPE diagram which accounts for about 50% of 
the empirical Acicib, Eq. (|3.6p, (cf. Table |TV|) . 

In pp scattering, the one-pion-exchange potential, V OPE , is given by, 
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V OPE (pp) = Ko , (3.8) 

while in T = 1 np scattering, we have, 

V OPE (n Pl T = 1) = -V n o + 2V n ± . (3.9) 

If the pion masses were all the same, these would be identical potentials. However, due to 
the mass splitting, the T = 1 np potential is weaker as compared to the pp one. This causes 
a difference between T = 1 pp and np that is known as CIB. For completeness, we also give 
the T = np OPE potential which is 

V OPE (np, T = 0) = -V w o - 2V n ± . (3.10) 

Due to the small mass of the pion, OPE is also a sizable contribution in all partial 
waves with L > 0; and due to the pion's relatively large mass splitting (3.4%), OPE creates 
relatively large charge-dependent effects in all partial waves (cf. Tables [V|and [VI] and Fig. [|). 
Therefore, all modern phase shift analyses |T^J4^] and all modern NN potentials 
include the CIB effect created by OPE. 

However, pion mass splitting creates further CIB effects through the diagrams of 2-7T 
exchange and other two-boson exchange diagrams that involve pions. The evaluation of 
this CIB contribution is very involved, but it has been accomplished in Ref. ||. The CIB 
effect from all the relevant two-boson exchanges (TBE) contributes about 1.3 fm to AaciB- 
Concerning phase shift differences, it is noticable up to D waves and can amount up to 50% 
of the OPE effect in some states (cf. Tables [V] and gjj). 

Another source of CIB is irreducible 717 exchange. Recently, these contributions have 
been evaluated in the framework of chiral perturbation theory by van Kolck et al. [[|(J . Based 
upon this work, we have calculated the impact of the 717 diagrams on the 1 5'o scattering 
length and on np phase shifts, (see column '^7' in Tables 0, 0, and 0). In L > states, 



the size of this contribution is typically the same as the CIB effect from TBE. 

In the 1 S'o state, the 717 contribution increases the discrepancy between theory and 
experiment (cf. Table |TV| ). As a matter of fact, about 25% of AaciB is not explained. For 
that reason, a quantitative fit of the empirical AaciB requires a small phenomenological 
contribution. The same is true for the difference between the empirical np and pp phase 
shifts in the 1 S , state (cf. Table [V|). 

For convenience, the major CIB effects on the strong NN force are plotted in Fig. In 
Fig. |5] the total CIB phase shift effect caused by the strong force is compared to the Coulomb 
effect on pp phase shifts (5 C denotes the phase shift in the presence of the Coulomb force, 
see Appendix A. 3 for precise definitions of 5 and 5 ). 

From the figures and tables it is evident that TBE and 717 create sizable CIB effects in 
states with L > 0. Therefore, we will include these two effects in our model. We note that 
conventional charge- dependent iViV models ignore these two contributions. 



IV. NUCLEON-NUCLEON SCATTERING 

We construct three NN interactions: a proton-proton (pp), a neutron-neutron (nn), and 
a neutron-proton (np) potential. The three potentials are not independent. They are all 
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based upon the model described in Sec. II and the differences between them are determined 
by CSB and CIB as discussed in Sec. III. Thus, when one of the three potentials is fixed, 
then the T = 1 parts of the other two potentials are also fixed due to CSB and CIB. 

We start with the pp potential since the pp data are the most accurate ones. Data fitting 
is done in three steps. In the first step, the pp potential is adjusted to reproduce closely 
the pp phase shifts of the Nijmegen multi-energy pp phase shift analysis f42|. This is to 
ensure that phase shifts are in the right ballpark. In the second step, the x 2 that results 
from applying the Nijmegen pp error matrix jPJ is minimized. The error matrix allows to 
calculate the x 2 in regard to the pp data in an approximate way requirings little computer 
time. Finally, in the third and crucial step, the pp potential parameters are fine-tuned by 
minimizing the exact x 2 that results from a direct comparison with all experimental pp 
data. During these calculations, it was revealed that the Nijmegen pp error matrix yields 
very accurate x 2 U P to 75 MeV. Therefore, in this final step, we used the error matrix up to 
75 MeV and direct x 2 calculations above this energy. 

The nn potential is constructed by starting from the pp potential, leaving out the 
Coulomb force, changing the nucleon masses, and fine-tuning the parameters such that 
the CSB differences listed in Tables |I| and |II| are reproduced. 

Concerning the np potential, we need to distinguish between the T = and T = 1 
states. In T = 1, we start from the pp potential, leave out the Coulomb force, change the 
nucleon masses, and replace the pp OPE potential by the one that applies to np. This 
produces a large part of CIB. The additional CIB due to TBE and discussed in Sec. Ill 
is incorporated by fine-tuning the parameters such that the total CIB phase shift differences 
as given in Table |VT] are reproduced. This fixes the np potential in the T = 1 states with 
L > 0. The 1 S'o np potential is adjusted such as to minimize the x 2 i n regard to the np 
data. The np T = potential is fixed by going through the entire three-step procedure: fit 
of Nijmegen T = phase shifts, minimizing the approximate x 2 obtained from the Nijmegen 
error matrix, and finally minimizing the exact x 2 that results from a direct comparison with 
all experimental np data. 

The resulting phase shifts for pp, nn, and np scattering in partial waves with J < 4 are 
given in Tables |VU| - [X|; pp phase shifts are plotted in Fig. |6| and np phase shifts are shown 
in Fig. |7|. For pp scattering, we show the phase shifts of the nuclear plus relativistic Coulomb 
interaction with respect to Coulomb wave functions; that is — in the notation of Ref . - 
we use Vc = ot' /r for the Coulomb potential and calculate the phase shifts Sq +n (= 5 C in 
our notation). We note that, for the calculation of observables (e. g., to obtain the x 2 i n 
regard to experimental data), we use electromagnetic phase shifts, as necessary, which we 
obtain by adding to the Coulomb phase shifts the effects from two-photon exchange, vacuum 



polarization, and magnetic moment interactions as calculated by the Nijmegen group [46,47 



This is important for 1 So below 30 MeV and negligible otherwise. For nn and np scattering, 
we show the phase shifts of the nuclear interaction with respect to Riccati-Bessel functions. 
All details of our phase shift calculations are given in Appendix A. 3 

The low-energy scattering parameters are shown in Table KB. For nn and np, the effective 
range expansion without any electromagnetic interaction is used. In the case of pp scattering, 
the quantities and are obtained by using the effective range expansion appropriate in 
the presence of the Coulomb force (see Appendix A. 4 for details). Note that the empirical 
values for a^ p and that we quote in Table £X1] were obtained by subtracting from the 
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corresponding electromagnetic values the effects due to two-photon exchange and vacuum 
polarization. Thus, the comparison between theory and experiment conducted in Table XI 
is adequate. 

For the comparison with the NN data, we consider three databases: 1992 database, 
after-1992 data, and 1999 database. The 1992 database is identical to the one used by the 
Nijmegen group for their phase shift analysis f49"]42" |. It consists of all NN data below 350 
MeV published between January 1955 and December 1992 that were not rejected in the 
Nijmegen data analysis (for details of the rejection criteria and a complete listing of the 
data references, see Refs. pfi||f9|J32[] ). The 1992 database contains 1787 pp data and 2514 np 
data. 

After 1992, there has been a fundamental breakthrough in the development of exper- 
imental methods for conducting hadron-hadron scattering experiments. In particular, the 
method of internal polarized gas targets applied in stored, cooled beams is now working 
perfectly in several hadron facilities, e. g., IUCF and COSY. Using this new technology, 
IUCF has produced a large number of pp spin correlation parameters of very high precision. 
In Table |X11| , we list the new IUCF data together with other pp data published between 
January 1993 and December 1999. Table [XII| lists all published after-1992 pp data below 
350 MeV except for one set, namely, 14 pp differential cross sections at 45 deg (lab.) be- 
tween 299.8 adn 406.8 keV by Dombrowski et al. [56]; according to the Nijmegen rejection 



criteria [|||, this set is to be discarded. The total number of (accepted) after-1992 pp data 
is 1145, which should be compared to the number of pp data in the 1992 base, namely, 1787. 
Thus, the pp database has increased by about 2/3 since 1992. The importance of the new 
pp data is further enhanced by the fact that they are of much higher quality than the old 
ones. 

Neutron-proton data published between January 1993 and December 1999 and included 
in our \ 2 calculations are listed in Table [XIII| . There are 544 such data, which is a small 
number as compared to the 2514 np data of the 1992 base. Note that Table |XJ11| is not a list 
of all np data published after 1992. Not listed are four measurements of np differential cross 
sections |67]-f7(J. We have examined these data and found in each case that they produced 
an improbably high y 2 when compared to current phase shift analyses fl4"2| , |4"5|| . Applying 
the Nijmegen rejection rule j|6],f|2|, the data of all four experiments are to be discarded. 
We follow this rule here, because we use the Nijmegen database for the pre-1993 period. 
When we add data to this base, then consistency requires that we apply the same selection 
criteria used for assembling the older part of the base. However, we like to stress that we do 
understand that any discarding of published data (i. e., data that have passed the refereeing 
process) is a highly questionable procedure. The problem of the np differential cross section 
data is an unresolved issue that deserves the full attention of all NN practitioners. Some 
aspects of the problem were recently discussed in Ref. |7l| . 



Finally, our 1999 database is the sum of the 1992 base and the after-1992 data and, thus, 
consists of the world NN data below 350 MeV that were published before the year of 2000 
(and not rejected). 

The xVdatum produced by the CD-Bonn potential in regard to the databases defined 
above are listed in Table [XI V] For the purpose of comparison, we also give the corresponding 
X 2 values for the Nijmegen phase shift analysis |^2| and the recent Argonne V\8 potential fl3~2| 



What stands out in Table |X1V| are the rather large values for the % 2 /datum generated by 
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the Nijmegen analysis and the Argonne potential for the the after-1992 pp data, which are 
essentially the new IUCF data. This fact is a clear indication that these new data provide a 
very critical test /constraint for any NN model. It further indicates that fitting the pre-1993 
pp data does not nessarily imply a good fit of those IUCF data. On the other hand, fitting 
the new IUCF data does imply a good fit of the pre-1993 data. The conclusion from these 
two facts is that the new IUCF data provide information that was not contained in the old 
database. Or, in other words, the pre-1993 data were insufficient and still left too much 
lattitude for pinning down NN models. One thing in particular that we noticed is that the 
3 Pi phase shifts above 100 MeV have to be lower than the values given in the Nijmegen 
analysis. 

The bottom line is that, for the 1999 data base (which contains 5990 pp and np data), 
the CD-Bonn potential yields a x 2 /datum of 1.02, while the Nijmegen analysis produces 
1.04 and the Argonne potential 1.21. We have also compared other recent NN potentials 
and NN analyses to the 1999 database and found in all case a x 2 /datum > 1.05. 

Thus we can conclude that the CD-Bonn potential fits the world NN data below 350 
MeV available in the year of 2000 better than any phase shift analysis and any other NN 
potential. 



V. THE DEUTERON 

The CD-Bonn potential has been fitted to the empirical value for the deuteron binding 
energy Bd = 2.224575 MeV [[72] using relativistic kinematics. Once this adjustment has 



been made, the other deuteron properties listed in in Table |XV| are predictions. For the 



asymptotic D/S state ratio, we find rj = 0.0256 — in accurate agreement with the empirical 
determination by Rodning and Knutson WM . The deuteron matter radius is predicted to be 



rd = 1.966 fm which agrees well with the value extracted from recent hydrogen-deuterium 



isotope shift measurements, r d = 1.971(6) fm []75]. Note that the deuteron effective range 
p d = p(—B d , —Bd) and the asymptotic S state A s are not directly observable quantities. 
Thus, 'empirical' values for p d and A s quoted in the literature are model dependent. There- 
fore, the perfect agreement between our predictions and the empirical values for p d and A s 
is of no fundamental significance. It only means that all models (including our own) are 
consistent with each other. 

More interesting is our prediction for the deuteron quadrupole moment Qd = 0.270 
fm 2 which is below the empirical value of 0.2859(3) fm 2 [[^,[7^]. Our calculation does not 



include relativistic and meson current corrections which according to Henning |77| contribute 
typically about 0.010 fm 2 for the Bonn OBE potentials. This would raise our theoretical 
value to Qd ~ 0.280 fm 2 , still 0.006 fm 2 below experiment. All recent NN potentials that 
use the 'small' irNN coupling constant g 2 /47r = 13.6 underpredict Qd by about the same 
amount. In Refs. |78|Ji7[| it was shown that Qd depends sensitively on g n and that a value 
<7 2 / 47T > 14.0 would solve the problem. However, a larger g n is inconsistent with the low- 
energy pp A y data (see Ref. [|I7|] for a detailed discussion of this issue). Thus, the accurate 
explanation of the deuteron quadrupole moment is an unresolved problem at this time. 

In Table IXVl we also give the deuteron D-state probability Pq. This quantity is not an 
observable, but it is of great theoretical interest. CD-Bonn predicts Pd = 4.85% while local 
potentials typically predict Pp « 5.7 %, which is clearly reflected in the deuteron D- waves, 
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Figs. P and |9[ The smaller value of CD-Bonn can be traced to the non-localities contained 
in the tensor force as discussed in Sec. II and demonstrated in Fig. f|. The CD-Bonn and 
the Nijmegen-I |3T[ potentials have nonlocal central forces which explains the soft behavior 



of their deuteron S'-waves at short distances that is particularly apparent in the plot of 
Fig. [| Numerical values of our deuteron waves and a convenient parametrization are given 
in Appendix D which also contains an account of how to conduct deuteron calculations in 
momentum space. 

VI. CONCLUSIONS 

We have constructed charge-dependent NN potentials, that fit the world proton-proton 
data below 350 MeV (2932 data) with a xVdatum of 1.01 and the corresponding neutron- 
proton data (3058 data) with % 2 /datum = 1.02. This reproduction of the NN data is more 
accurate than by any other known NN potential or phase-shift analysis. A particular chal- 
lenge are the pp spin correlation parameters that were recently measured at the IUCF Cooler 
Ring with very high precision (1126 data below 350 MeV). Our pp potential reproduces these 
data with x 2 /datum = 1.03, while the high-quality Nijmegen analysis |^2] and the Argonne 
Vis potential produce x 2 /datum of 1.24 and 1.74, respectively, for these data. 

The charge-dependence of the present potential (that has been dubbed 'CD-Bonn') 
is based upon the predictions by the Bonn Full Model for charge-symmetry and charge- 
independence breaking in all partial waves with J < 4. Thus, our model includes con- 
siderably more charge- dependence than other recently developed charge-dependent poten- 
tials |3"T,|3"2lj. For example, the Nijmegen potentials |JT] include essentially only charge- 



dependence due to OPE which produces CIB, but no CSB. Thus, the Nijmegen group does 
not offer any genuine neutron-neutron potentials. To have distinct pp and nn potentials is 
important for addressing several interesting issues in nuclear physics, like the 3 H- 3 He bind- 
ing energy difference for which the CD-Bonn potential predicts 60 keV in agreement with 
empirical estimates. Another issue is the Nolen-Schiffer anomaly [If. Some potentials that 
include CSB focus on the 1 S'o state only, since this is where the most reliable empirical 
information is. However, this is not good enough. In Ref. P it has been shown that CSB 
in states with J > is crucial for the explanation of the Nolen-Schiffer anomaly. 

The CD-Bonn potential is represented in terms of the covariant Feynman amplitudes for 
one-boson exchange which are nonlocal. Therefore, the off-shell behavior of the CD-Bonn 
potential differs in a characteristic and well-founded way from the one of commonly used 
local potentials. 

The simplest system in which off-shell differences between NN potentials can be in- 



vestigated is the deuteron (see Ref. |79| for a thorough study of this issue). Our plots of 
the deuteron wave functions, Figs. H and make this point very clear. Empirical tests 
of deuteron wave functions can be conducted via the structure functions A(Q 2 ), B{Q 2 ), 
and the tensor polarization in elastic electron-deuteron scattering T 2 o(Q 2 ) or, alternatively, 
via the three deuteron form factors Gc(Q 2 ), Gq(Q 2 ), and Gm(Q 2 ), for which the deuteron 
wave functions are crucial input. Using the deuteron wave functions derived from the Bonn 
model, Arenhovel and coworkers |SIJ find a good agreement between theory and experiment 
for A(Q 2 ), B(Q 2 ), and T 20 (Q 2 ) up to Q 2 = 30 fm -2 . Very recently, the tensor polarization 
^2o(Q 2 ) nas been measured up to Q 2 = 45 fm -2 at the Jefferson Laboratory \WL\. The best 
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reproduction of these new high-precision data is provided by two calculations that are based 



upon the Bonn deuteron wave functions [p2| , |33 



Another way in which the off-shell behavior of our potential shows up is by yielding larger 
binding energies in microscopic calculations of nuclear few- and many-body systems |84 



where underbinding is a persistent problem. To demonstrate this, we have computed the 
binding energy of the triton in a 34-channel, charge- dependent Faddeev calculation. The 
prediction by the CD-Bonn potential is 8.00 MeV. Local potentials typically predict 7.62 
MeV |85|,|86| and the experimental value is 8.48 MeV. Thus, the nonlocality of the CD-Bonn 
potential explains almost 50% of the gap that persists between the predictions by local 
potentials and experiment. Similar results are obtained for the a particle JSBJSTjj. Concerning 



the small difference that is left between the CD-Bonn predictions and experiment, two 
comments are in place. First, besides the relativistic, nonlocal effects that can be absorbed 
into the two-body potential concept, there are further relativistic corrections that come from 
a relativistic treatment of the three-body system. This increases the triton binding energy 
by 0.2-0.3 MeV |88|-pO|,P|. Second, notice that the present nonlocal potential includes only 



the nonlocalities that come from meson-exchange. However, the composite structure (quark 



substructure) of hadrons should provide additional nonlocalities |91[ which may be even 
larger. It is a challenging topic for future research to derive these additional nonlocalities, 
and test their impact on nuclear structure predictions. 

The trend of the nonlocal Bonn potential to increase binding energies has also a very 
favorable impact on predictions for nuclear matter |7|p2]| and the structure of finite nu- 
clei m-g. 



Due to the very accurate fit of even the latest high-precision NN data; due to the 
comprehensive and sophisticated charge- dependence incorporated in the model; and due 
to the well-founded off-shell behavior, the CD-Bonn potential [^3J represents a promising 
starting point for exact few-body calculations and microscopic nuclear many-body theory. 
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APPENDIX A: TWO-NUCLEON SCATTERING IN MOMENTUM SPACE 

1. Scattering equation 



Two-nucleon scattering is described covariantly by the Bethe-Salpeter (BS) equation |2T 
which reads in operator notation 

T = V + VQT (Al) 

with T the invariant amplitude for the two-nucleon scattering process, V the sum of all 
connected two-particle irreducible diagrams, and Q the relativistic two-nucleon propaga- 
tor. Since this four- dimensional integral equation is very difficult to solve, so-called three- 
dimensional reductions have been proposed, which are more amenable to numerical solution. 



15 



Furthermore, it has been shown by Gross [96[] that the full BS equation in ladder approx 



imation (that is, the kernel V is restricted to the exchange of single particles as, e. g., in 
the OBE model) does not have the correct one-body limit (i. e., when one of the particles 
becomes very massive) while a large family of three-dimensional quasi-potential equations 
does. These approximations to the BS equation are also covariant and satisfy relativistic 
elastic unitarity. Three-dimensional reductions are typically derived by replacing Eq. ( |A1| ) 
with two coupled equations P7| : 



T = W + WgT 



(A2) 



and 



yy = V + V(G - g) W 



(A3) 



where g is a covariant three-dimensional propagator with the same elastic unitarity cut as 
Q in the physical region. In general, the second term on the r.h.s. of Eq. ( |A5| ) is dropped to 
obtain a true simplification of the problem. 

More explicitly, the BS equation for an arbitrary frame reads [ITS] 



T(q'; q\P) = V(q'; q\P) + J d 4 k V(q'; k\P) Q(k\P) T(k; q\P) 



with 



g(k\p) 



i 



i 



a p+ 



% 

2^ 



M + ie)« (| P- ^-M + ie 
fi + M 1 (1) [ \ P- 



1(2) 



k) 2 -M 2 + ie 



[\P — k) 2 — M 2 + ie 



(2) 



(A4) 

(A5) 
(A6) 



where q, k, and q' are the initial, intermediate, and final relative four-momenta, respectively, 
and P = (Pq, P) is the total four-momentum. For example, in the initial state we have: 
Q — \{Vi ~ P2), P = Pi +P2, and pi/ 2 = \P ± q with p 1 and p 2 the individual four-momenta 
of particle 1 and 2. In the center-of-mass (cm.) frame, we will have P = (-y/s, 0) with yfs 
the total energy. For all four-momenta, our notation is k = (fco,k); \jt = 7^^. M denotes 
the nucleon mass. The superscripts in Eq. ( |A6| ) refer to particle (1) and (2). At this stage, 
T, V, and Q are operators in spinor space, i. e. they are 16 x 16 matrices which, when 
sandwiched between Dirac spinors, yield the corresponding matrix elements. 

It is common to the derivation of all three-dimensional reductions that the time compo- 
nent of the relative momentum is fixed in some covariant way, so that it no longer appears 
as an independent variable in the propagator. 

Following Blankenbecler and Sugar (BbS) |[20| , one possible choice for g is (stated in 
manifestly covariant form for an arbitrary frame) 



gBbs{k, s) 



ds' 



Um 2 s' — s — ie 
x5 (+) [(iP' -k) 2 -M 



£(+)[ (lp' + k) 2 — M 2 



x 



\ P'+ ^ + M]«[| P>- ji + M 



(2) 



(A7) 
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with indicating that only the positive energy root of the argument of the 5-function is 
to be included; P 2 = s and P' = ^P- By construction, the propagator g BbS has the same 
imaginary part as Q and, therefore, preserves the unitarity relation satisfied by T. In the 
cm. frame, integration yields 

9Bbs(k, s) = 6(ko) g B bs(K s) (A8) 

with 

„ N M 2 Aj } (k) Ai 2) (-k) /A , 

g BbS (k, s) = — +i (A9) 



AS?(k) = f ^-:; k + M | (A10) 



where 

(i) 

IS. -f XVI \ 

2M / 

= X;i«(k,A i )>(ii(k ) A i )| (All) 

represents the positive-energy projection operator for nucleon i {% = 1 or 2) with ii(k) a 
positive-energy Dirac spinor of momentum k; u = m^7°. A« denotes the helicity of the 
respective nucleon, and £7*. = \/ M 2 + k 2 with M the nucleon mass. The projection operators 
imply that virtual anti-nucleon contributions are suppressed. 

Using the approximation W ~ V [cf. Eq. (|A3|) 1, we obtain the explicit form of Eq. ( |A2|) 
by simply replacing Q by g B bs i n Eq. (|A4j). This yields in the cm. frame 



T(0, q'; 0, q| y/s) = V(0, q'; 0, q) + J d 3 k V(0, q'; 0, k) g BbS (k, s) T(0, k; 0, q\y/I). (A12) 

Note that four-momentum is conserved at each vertex, and that in the initial state the 
nucleons are on their mass-shell, therefore q = (0, q). The total cm. energy is 



4~s = 2E q = 2^/M 2 + q 2 . (A13) 
With this we obtain, simplifying our notation, 

T(q', q) = V(q', q) + / Sk V(q', k) ^ ^ ^ ^ ^ T(k, q). (A14) 

Taking matrix elements between positive-energy spinors yields an equation for the invariant 
scattering amplitude 

r M 2 1 

r(q', q) = V(q, q) + / d 3 * V(q', k) — 2 f (k, q), (A15) 

J hib Q — k z + ze 



where helicity and isospin indices are suppressed. 
Defining 



T(q / ,q) = ^T(q',q)^ (A16) 
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and 



V(q',q) 



~rr ^(q ! q 




(A17) 



which has become known as "minimal relativity" [98], we can rewrite Eq. ( ]A15| ) as 

T(q', q) = V(q', q) + / d 3 * V(q', k) 2 f 2 . T(k, q) 

J q^ — k^ + ze 



(A18) 



which has the form of the familiar Lippmann-Schwinger equation. The quantity T has the 
usual (nonrelativistic) relation to phase shifts and AW observables. Thus, the AW potential 
V defined in Eq. (|A17|) and used in the above Lippmann-Schwinger equation can be applied 
in the deuteron and in conventional nuclear structure physics in the same way as any other 
(nonrelativistic) potential. This is the great virtue of the (relativistic) BbS equation. 



2. R-matrix and partial wave decomposition 

In solving the scattering equation, it is more convenient to deal with real quantities. We 
shall therefore introduce the real i?-matrix (better known as 'A'-matrix') defined by 

T = R- inR5(E - H Q )T (A19) 

The equation for the real i?-matrix corresponding to the complex T-matrix of Eq. ( |A18| ) is 

i?(q', q) = V(q[, q)+vj d 3 k V(q[, k) J2(k, q) (A20) 

where V denotes the principal value. 

Now, we need to also include the spin of the nucleons. Relativistic scattering of particles 



with spin is treated most conveniently in the helicity formalism [99]. Therefore, we will use 
a helicity state basis in our further formal developments. Our presentation will be relatively 
brief; a more detailed derivation is given in Appendix C of Ref. U which is based upon 
Refs. PlJlOOj. 



The helicity Aj of particle i (with % = 1 or 2) is the eigenvalue of the helicity operator 
|<Tj • which is ±§. 

Using helicity states, the i?-matrix equation reads, after partial wave decomposition, 

(KK\ rJ (l' i <?)|AiA 2 ) = (X' 1 X' 2 \V J (q' } q)\\ 1 \ 2 ) 

+ Y,V dkk 2 ^-—{W 2 \V J {q',k)\h x h 2 ) 



hi ^2 



x (/i 1 /i 2 | J R J (A;,g)|AiA 2 ) (A21) 

where J denotes the total angular momentum of the two nucleons. Here we are changing 
our notation for momenta: in the above equation and throughout the rest of Appendix 
A, momenta denoted by non-bold letters are the magnitude of three-momenta, e. g. q = 
|q|, k = |k|, etc.; hi and h 2 are the helicities in intermediate states for nucleon 1 and 2, 
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respectively. Equation ( |A21|) is a system of coupled integral equations which needs to be 
solved to obtain the desired matrix elements of R J . 

Ignoring anti-particles, there are 4 x 4 = 16 helicity amplitudes for R J . However, time- 
reversal invariance, parity conservation, and the fact that we are dealing with two identical 
fermions imply that only six amplitudes are independent. For these six amplitudes, we 
choose the following set: 



R( (</, q) = {+ + 


\R J (q',q)\ 


+ +) 


Ri(q',q) = (+ + 


\R J (q',q)\ 




R J 3 (q',q) = (+- 


\R\q\q) 


+ -) 


R J Aq',q) = (+- 


\R J (q',q) 


- +) 


R J 5 (q',q) = (+ + 


\R J (q',q)\ 


+ -) 


Ri(q,q) = (+- 


\R\q\q) 


+ +) 


Notice that 






Rt(q',q) = 







(A22) 



(A23) 

We have now six coupled equations. To partially decouple this system, it is usefull to 
introduce the following linear combinations of helicity amplitudes: 



(A24) 



°R J 


= n x — n 2 


l R J 


= K 3 ~ K A 


12 R J 


= R"[ + i?2 


M R J 


= -Rg + -R4 


55 R J 


= 2Rl 


m R J 


= 2Ri 



We also introduce corresponding definitions for V J . Using these definitions, Eq. (|A21|) 
decouples into the following three sub-systems of integral equations: 
Spin singlet 



M 



°Rj(q>, q ) = °V J (q', q)+V J dk k 2 -f\^ °V J (q 1 , k) °R J (k, q) 
Uncoupled spin triplet 

/■EX) M 

l R J {q', q) = V J (g', q)+V dk k 2 — — V J (g', k) l R J (k, q) 

Jo q z — k z 

Coupled triplet states 



(A25) 



(A26) 



roc 

12 R J {q', q) = 12 V J (q', q)+V dk k 2 

Jo 

+ 55 V J (q',k) m R J (k,q)} 

roo 

R J (q\ q ) = u V J (q', q) + V / dk k 2 

Jo 



M 



34 D J 



q 2 - k< 
M 



i 12 V J (q',k) 12 R J (k,q) 



[ M V J {q',k) ' M R J (k,q) 



34 ryJi 



k 2 
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+ m V J {q',k) 55 R J {k,q)} 



roc TUT 

55 R J (q', q) = 55 V J (q\ q)+P dk k 2 --—\ 12 V J (q', k) 55 R J (k, q) 

JO Q — K 



+ 55 V J (q\k) 3 *R J (k : q)} 



r°° M 

66 R J (q>, q) = ™V J (q', q)+V dk k 2 — -[ U V J (q', k) m R J (k, q) 

Jo q z — k z 



+ m V J (q',k) l2 R J (k,q)}. 



(A27) 



More common in nuclear physics is the representation of two-nucleon states in terms of 
an \LSJM) basis, where S denotes the total spin, L the total orbital angular momentum, 
and J the total angular momentum with projection M. In this basis, we will denote the R 
matrix elements by Rjf L = (L'SJM\R\LSJM). These are obtained from the helicity state 
matrix elements by the following unitary transformation: 
Spin singlet 



Uncoupled spin triplet 



Coupled triplet states 



pJO 



°R J 



l R J 



(A28) 



(A29) 



p71 

"-7-1,7- 




1 




1 = 


2J + 


1 


p71 

"-74-1,74 




1 




-1 = 


2J + 


1 


p71 

•"■7-1,74 




1 




-1 = 


2J + 


1 


p71 

"■74-1,7- 




1 




1 = 


2J + 


1 



J 12 R J + (J + 1) M R J + ^J(J + 1)( 55 R J + m R J ) 
(J + 1) 12 R J + J U R J - v /j(J+l)( 55 R J + m R J ) 
J(J+ 1)( 12 R J - U R J ) - J 55 R J + ( J + 1) m R J 



J {J + 1)( 12 R J - U R J ) + (J + 1) bb R J - J m R J ) 



55 pJ 



(A30) 



Similar notation and transformations apply to V. 

One way to proceed is to solve the system of equations ( A27 ) and then apply the trans- 
formation ([A30|) . Alternatively, one may apply the transformation ( |A30| ) directly in ( |A27| ) 
to obtain the system of four coupled integral equations in LSJ representation, 



R J + 1 + (q',q) = V^ + (q',q)+V 



dk k< 



M 



[V'Hq',k) Rt\{Kq) 



71 



q 2 -k 2i ' +H 



roc TUT 

R n -(q', q) = V n _(q', q)+V dkk 2 — -[ Vfitf, k) R J ±(k, q) 

Jo q 2 — k 



+ V^(q',k)Ril(k,q)} 



71 



R J + 1 _(q',q)=V^(q',q)+V f 

J 



"dk*-" 



Vf + (q',k)R J + l Jk,q) 



71 



k 2 



+ V£{q',k)R J _\{k,q)\ 



71 
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roo A/T 

R-+tf, Q) = V^ + (q', q)+-PJ o dk e - T -^[ V n _{q', k) R J _\(k } q) 



+ V n + {q\k)R J + \{k,q)] 



M 



(A31) 



where we used the abbreviations R+ + = Rj\u+x, R- 1 - = Rj-u-ii R J +- = 
Rj+i,j-i, Rl\ = Rj-i,j+i, and similarly for V. 

The above integral equations can be solved numerically by the matrix inversion 
method [|101|| . The method is explained in detail in Ref ||102|| where also a computer code is 
provided. 

Each two-nucleon state carries a well-defined total isospin T (which is either or 1) that 
is fixed by 

(_1)^+s+t = _ x (A32) 



3. Phase shifts 

Phase shifts are determined from the on-energy-shell i?-matrix through: 
Spin singlet 

tan °5 J (T lab ) = ~qM °R J (q, q) (A33) 

Uncoupled spin triplet 

tan 1 5 J (Ti ab ) = ~qM 1 it! J (g, q) (A34) 

For the coupled states, a unitary transformation is needed to diagonalize the two-by-two 
coupled .R-matrix. This requires an additional parameter, known as the 'mixing parameter' 
ej. Using the convention introduced by Blatt and Biedenharn [|103j| , the eigenphases for the 
coupled channels are, in terms of the on-shell i?-matrix, 

tan 8 J (T lah ) 



tan2e j(Ti ab ) 

Here, all i?-matrix elements carry the arguments (q, q) where q denotes the cm. on-energy- 
shell momentum. Based upon correct (relativistic) kinematical considerations, the momen- 
tum q and the nucleon mass M to be used in the above formulae are determined to be: 
Proton-proton scattering: 

q 2 = ^M p T lab , (A36) 
M = M p . (A37) 



7T 



--qM 



+ Rj+i,j+i ^ 



R j j- 



i,j-i 



-"-J+l.J+l 



cos2ej 



2Rj+i,j-i 



pj 

n j-i,j-i 



n j+i,j+i 



(A35) 
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Neutron-neutron scattering: 



q 

M 



-M n T lab 

Mr,.. 



Neutron-proton scattering: 



Q 

M 



M*T lab (T lab + 2M n ) 
(M p + M n f + 2T lab M p 
2M p M n 



M p + M n 



(A38) 
(A39) 

(A40) 
(A41) 



In the above, M p denotes the proton mass, M n the neutron mass (see Table | for their 
accurate numerical values) and Ti ab is the kinetic energy of the incident nucleon in the 
laboratory system. 



An alternative convention for the phase parameters has been used by Stapp et al. ||104 
known as 'bar' phase shifts. These are related to the Blatt-Biedenharn parameters by 



sin(^ 
sin {8 J _ 



5 J + + 5 J _ 

si: 



5 J + + St 



tan2ej/tan2ej 
sin 2ej/sin 2ej 



(A42) 



In this paper, all phase shifts shown in tables or figures are in the 'bar' convention, even 
though we omit the bar in our notation. 

The above formulae apply to the calculation of phase shifts when only the short-range 
nuclear force is taken into account (and no electromagnetic interaction). This is, in general, 
appropriate for nn and np scattering. We also note that the above momentum space method 
is exactly equivalent to calculations conducted in r-space where the radial Schrodinger equa- 
tion, 



f_ 2 

dr 2 



L{L + 1) 



-MV 



XL (r; q) = 



(A43) 



is solved for the radial wave function Xz(r;g) which is then matched to the appropriate 
asymptotic form of the wave function to obtain the phase shift. When no long-range potential 



is involved, the asymptotic wave functions are Riccati-Bessel functions |j!05| . 

In pp scattering, the long-range Coulomb potential must be taken into account. The 
asymptotic form of the wave function then is 



Xi(r;q) oc F L (r)',qr) + tan 6 L G L (r)',qr) 



(A44) 



with Fl and Gl the regular and irregular Coulomb functions |1U5 |. By 5 C we denote the 
phase shift of the nuclear plus Coulomb interaction with respect to Coulomb wave functions; 
that is, in the notation of Ref. f46fl , 5 C = S^ +N . The parameter rj' is the "relativistic" rj 
defined by TOFFW 



Vlab 



Mp 

2q 



-ft 



(A45) 
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with 

E 2 + a 2 

a' = a « „ , (A46) 

MpE q 1 > 

and a = 1/137.035989 ||. The total potential V that appears in Eq. ( |A43| ) is now the sum 
of the nuclear potential Vn and the Coulomb potential Vc] i. e., 

V = V N + V C , (A47) 

where we use the "relativistic" Coulomb potential ||107 



V c = - . (A48) 
r 

Since we conduct our calculations in momentum space, we do not solve Eq. ( |A43| ) and, thus, 
do not have a numerical x( r ; <?) available that can be matched directly to the asymptotic 
form Eq. ( |A44|) . However, there are ways to perform this matching within the framework 
of momentum space calculations. We follow here the method proposed by Vincent and 



Phatak ||108|| in which the potential is divided into a short-range part V$ and a long-range 

V = V S + V L (A49) 



part Vl] i- e., 



with 

V s =(y N + Vc)9(R-r), (A50) 
V L = V c 6(r - R) , (A51) 

where R is to be chosen such that the short-range nuclear potential has vanished for r > R 
(R ps 10 fm is an appropriate choice); and 9 is the usual Heaviside step function. First, 
one calculates the phase shift (denoted by <5f ) that is produced by Vg alone. Notice that 
Vs is of range R and consists of the nuclear potential plus the Coulomb potential cut off at 
r = R. There is no problem in performing numerically the Bessel transformation of a cutoff 
Coulomb potential to produce the momentum space version of this potential for the various 
partial waves. Since Vs is of finite range, the momentum space formalism can be used to 
calculate 5f . The asymptotic wave function associated with Vs and <5f is 

X S L (r; q) oc F L ( V ' = 0, gr) + tan5 s L G L (r]' = 0, qr) (A52) 

which should match smoothly the asymptotic function Eq. (|A44j) at r = R. Note that 
F L (rj' = 0, gr) and Gl(t)' = 0, gr) are equal to Riccati-Bessel functions. Matching the 
logarithmic derivatives yields the desired formula for the phase shift 

t»n/P - +tBn6 s L [G L (0),F L (rf)] 

L [G L ( V '),F L (0)]+Un5i[G L ( V '),G L (0)} ' [ ^ 6) 

where the square brackets denote the Wronskian 
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(A54) 

r=R 

and Fl(0) = Fl(t)' = 0, qr), Fl^t]') = Fl{t]' ■, qr); similarly for Gl- 

All pp phase shifts shown in this paper are Coulomb phase shifts, 5 C , as defined and 
calculated above. However, we like to stress that, for the calculation of observables (e. g., 
to obtain the x 2 i n regard to experimental data), we use electromagnetic phase shifts, as 
necessary, which we obtain by adding to the Coulomb phase shifts the effects from two- 
photon exchange, vacuum polarization, and magnetic moment interactions as calculated 
by the Nijmegen group |46|j47[1 . This is important for 1 S'o below 30 MeV and negligible 
otherwise. 



[F L (0),F L ( V >)) = 



dF L (j) 
dr 



- F L (rf 



dF L (Q) 
dr 



4. Effective range expansion 

For low-energy 5-wave scattering, q cot 6 can be expanded as a function of q 

"' q cot 5 w — + \rq 2 + C(g 4 ) (A55) 



tan 6 ' a 2 

where a is called the scattering length and r the effective range (for which, in some parts of 
this paper, we also use the notation a N and r ). This is appropriate for nn and np. 

In the case of pp scattering, where the Coulomb potential is involved, a more sophisticated 
effective range expansion must be applied [pfl , 



C 2 ( V ')qcot(6g) + 2q V 'h( V ') = -± + \r c pv q 2 + 0(q A ) , (A56) 

a pp L 

where Sp p denotes the 1 5'o pp phase shift with respect to Coulomb functions and Cq and h 
are the standard functions, 

CoW) = ^J^y , (A57) 
h(r/) = - ln(i/) + Re[4>(l + iri')\ (A58) 

OO 

= - ln(r/) - 7 + r, 12 ^ [n(n 2 + r/' 2 )]" 1 , (A59) 

n=l 

where ip denotes the digamma function and 7 = 0.5772156649 . . .. 

This formalism takes care of the Coulomb force. However, the full electromagnetic 
interaction between two protons has contributions beyond Coulomb, e. g., from two-photon 
exchange and vacuum polarization. To include the full electromagnetic interaction into 
the effective range expansion is very involved. Therefore, the empirical values for the pp 
effective range parameters (which naturally involve the full electromagnetic interaction) have 
been corrected (in a fairly model-independent way) for the electromagnetic effects beyond 
Coulomb ^J46|]. This procedure yields 'empirical' values for a pp and r pp which is what 



we quote in Table [X5 under 'Experiment'. The existence of empirical values of this kind 
makes the comparison between theory and experiment much easier; and it justifies that a 
theoretican calculates predictions for just a pp and r pp using the simple formalism outlined 
above. 
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APPENDIX B: ONE-BOSON EXCHANGE POTENTIAL 



1. OBE amplitudes 

The Lagrangians Eqs. ( |2.1|) -( [2.5D imply the following OBE amplitudes which we state 
here in terms of i times the Feynman amplitude: 

(q'AiA^|K|qA 1 A 2 ) 

<J; u(q', A;)z7 5 n(q, Ai)fi(-q', A 2 )i 7 5 «(-q A a )/[(q' - q) 2 + m*] , (Bl) 



(2tt / 

(q'A / 1 A / 2 |V;|qA 1 A 2 ) 

9l 



«(q', Ai)«(q, A 1 )«(-q / , A' 2 )«(-q, A 2 )/[(q' - q) 2 + m 2 ] , (B2) 



(2tt) 3 
(q'Ai^l^lqAiAa) 

:{fi(q , ) A' 1 ) 7 ^(q, Ai)}{u(-q / , A 2 ) 7 ^(-q, A 2 )}/[(q' - q) 2 + m 2 ] , (B3) 



(2vr) ; 



(q'AiA 2 |V;|qA 1 A 2 ) 



T2 -{^P M (q' ) A' 1 )7 M M(q, Ai) + -^-M(q, A'Jo^g' - g)"u(q> Ai)} 



(2tt) 3 ^ " v " 17 2M p 

2M„ 



x{g p u(-q', A 2 ) T ^(-q A 2 ) - ^(-q', A 2 K"i(g' - g),n(-q A 2 )} 



/[(q'-q) 2 + m 2 ] 

w(q', Ai)[(?' + g)p + (£' - £)(<7„o - 7„7o)Mq, A x )} 



2M P 



x{(g p + f p )u(-q', A' 2 )7"M(-q, A 2 ) 



2M p 



«(-q', A 2 )[(g' + gf + (£' - E)(g»° - T V)]«(-q, A 2 )} 



/[(q'-q) 2 + m 2 ], (B4) 

where for the pion we have suppressed isospin factors and charge-dependence which will be 
included later. Working in the two-nucleon cm. frame, the momenta of the two incoming 



(outgoing) nucleons are q and — q (q' and — q'). E = \JM 2 + q 2 , E' = J M 2 + q' 2 , and M is 
the nucleon mass. Using the BbS equation |20| , the four-momentum transfer between the two 
nucleons is (g'— g) M = (0, q' — q). The Gordon identity |H| has been applied in the evaluation 
of the tensor coupling of the p; (q 1 + g) M = (E' + E, q' + q) and (q' + q) 1 * = (E' + E, — q' — q) . 
The propagator for vector bosons is 

-9^ + {q' -?W -q)v/n% ^ 



(q' - q) 2 — m l 
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where we drop the {q' — q)^,{q ! — g)^-term which vanishes on-shell, anyhow, since the nucleon 
current is conserved. The off-shell effect of this term was examined in Ref. [ IPS ] and found 
to be unimportant. 

The Dirac spinors in helicity representation are given by 



u(q, Ai) 

u(~q, A 2 ) 
which are normalized such that 



l E + M ( 1 

2Ai]q] 
E+M 



2M 



! E + M 
2M 



1 

2A 2 |q| 
E+M 



|Ai), 
|A 2 >, 



(B6) 
(B7) 



u(q, A)u(q, A) = 1. 



(B8) 



with u = v)^° . 

At each meson-nucleon vertex, a form factor is applied which has the analytical form 



^[(q - q) 



A 2 -m 2 



a 2 + (q' - q) 



(B9) 



with m a the mass of the meson involved and A a the so-called cutoff mass. Thus, to obtain 
the final OBE potential V, the amplitudes Eqs. (|BT])-(|B~3|) are to be multiplied by J-^ and 
certain square-root factors [cf. Eq. (|2.7|) 1. 



2. Partial wave decomposition 



The potential is decomposed into partial waves according to 

(A' 1 A' 2 |F J (g',g)|A 1 A 2 ) = 2vr £ d{<x*0) d^_ M -%(WKK\V\&iXi) 

(BIO) 

where 6 is the angle between q and q' and d^ nm ,(9) are the conventional reduced rotation 
matrices which can be expressed in terms of Legendre polynominals Pj (cos 9). The following 
types of integrals occur: 



r(o) 



r(2) 



dt 



dt 



Pj(t) 



q) 2 + rn? a 



Qj(Zg) 

q'q 



tPj(t) 



Q 



q' 



q) 2 + m? a 



q'q 



J + l 



iJtPj{t) + Pj-i(t) Q 



dt 



(q' - q) 2 + m2 



q'q 



(Bll) 
(B12) 
(B13) 
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r(4) 



r(5) 



r(6) 



J 



+ 1 tPj(t) - Pj^(t) _Qf(z*) 



J + lJ-i " (q'-q) 2 + m 2 
t 2 Pj(t) Q (4) 



dt 



1 (^q'-q) 2 + m 2 
1 



-i ^ JfPjjfy+tP^) 
J + U-i (q' - q) 2 + m 2 



J 



J+lJ-i 



+i dt t 2 Pj(t)-tP J _ 1 (t) 



(q' - q) 2 + m 2 



q'q 
q'q 



(B14) 
(B15) 
(B16) 
(B17) 



with t = cos 9 and z Q = (q + q + rn a )/2q'q where our notation for momenta is q' = |q'|, 
and q = |q| which we will use throughout the remainder of the appendices. 

The Qj(z) are the Legendre functions of the second kind ||105|| ; e. g., Qo(z) = §ln[(z + 
l)/(z — 1)]. The combinations needed above are defined by: 



Q 



Q'f 



zQj - 5 jo , 
j^iJzQj + Qj-i) 



J 



J+l 



zQf - -Sji 



[zQj — Qj-\, 
1 



zQT - g*n 



zQf 



(B18) 
(B19) 

(B20) 

(B21) 

(B22) 

(B23) 
(B24) 



The integrals Eqs. ( [B11| )-( [BT7| ) can be evaluated either numerically or analytically by us- 
ing the Legendre functions of the second kind. The latter method is better if the correct 
threshold behavior of V J (q', q) for q', q — > is important. 

The above expressions still ignore the cutoff which is included by replacing 



1 



^[(q'-q) 



(q' 



q) 2 + m 2 



(q' 



q) 2 + m 2 



(B25) 



in Eqs. ( [B11| )-( |B~T7|) . If the Legendre functions of the second kind are used, then the product 
of propagator and cutoff must be decomposed according to 



^[(q'-q) s 



(q - qj 



mi 



'A 2 



mi. 



(q - q) 2 + m l 



A 2 



a 2 



A 2 aJ (q'-q) 2 + A 



+ 



Aai-m 2 



A 2 



a 2 



A 2 



al, 



(q - q 



2 + A 2 



(B26) 



q2 
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where A a i/ 2 = A a ± e with e — ► 0; i. e., e -C A a , e. g., e ~ 1 MeV. To give an example, if 
with cutoff is given by 

J( Q) [« dt Pj(t)J*[W-tf] (B27) 
J ~ J-i (q' - q) 2 + ml (B27) 
Qj(m Q ) / A 2 a2 -m 2 \ Qj(A al ) ( A 2 al - m\ \ Qj(A a2 ) 

and similarly for the other if. Notice that the if are functions of q', q, m a , and A a even 
though our notation does not indicate this. 



3. Final potential expressions 



Here, we will present the final potential expressions in partial wave decomposition. More 
details concerning their derivation can be found in Appendix E of Ref. |IJ. First, we state 
the potentials in terms of the combinations of helicity states defined in Eq. ( |A24j) . 



One-pion-exchange: 



with 



and 



OyJ 

' TV 

lyj 

TV 

12yJ 
34yJ 

' TV 

55yJ 

" 7T 

my J 



r(i)> 



a (fP if + *?> / 
a (-f<°> if - f« if) 

a if - Fi 0) if) 

CFflf 
-CFflf 



9l 



1 M M 
4vr 2ttM2 V^V^ 



= E'E - M 2 
= -q'q 
Ff = -M(E'-E). 



(B29) 



(B30) 



(B31) 



One-sigma-exchange : 



°K J = C a (F^lf+Fpif) 
% J = C„(FWlf+FPlf) 

12 V a J = C a (FPlf+Fflf) (B32) 
U Vj = CAF^lf+Fflf) 
5 V/ = C a Ff if 
m V a J = C a Fflf 
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with 



C„ 



9l 



<M M 



and 



4tt 2ttM 2 V E' V E 



Ff = -(E'E + M 2 ) 

= q'q 
Fl 2) = M(E' + E). 



One-omega-exchange: 



with 



°Vj = C w (2E'E - M 2 ) if 
l Vj^C u (E'Elf + q'qlf) 

12 V^C UJ (2q'qlf + M 2 lf) 

U Vj = C^ q ' q lf+E'Elf) 

55 <vj = -C u ME if 

<*Vj = -C u ME' if 



g 2 1 M 
u 4?r ttM 2 V E> 



The one-rho-exchange potential is the sum of three terms, 

v p = v vv + v tt + v vt . 



Vector-vector coupling 



with 



°V v l = C vv (2E'E - M 2 ) if 
% J V = C VV (E'E If + q'q If) 
l2 V v { = C m (2q'qlf+M 2 lf) 
U Kl = C vv (q'q if + E'E if) 

55 V v i = -Cw ME if 
™Vl = -C vv ME' if 



vv 



a 



vv 



91 t 1 -t 2 M M 

4tt ttM 2 V E'\ E 
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Tensor-tensor coupling 



°v t i = Cu{(q ,2 + q 2 )(3E'E + M 2 )if 

+ [q' 2 + q 2 - 2(3E'E + M 2 )]q'q if - 2q l2 q 2 if} 
V t f = C u {[Aq' 2 q 2 + {q' 2 + q 2 )(E' E - M 2 )] if 
+2(E'E + M 2 )q'q if 

-(q' 2 + q 2 + AE'E)q'q if - 2q' 2 q 2 if} (B40) 
i2 V t ! = C tt {[4M 2 -3(q' 2 + q 2 )}q'qlf 

+ [Qq' 2 q 2 - (q' 2 + q 2 )(E'E + 3M 2 )] if + 2(E'E + M 2 )q'q if} 
U V t } = C tt {-(q' 2 + q 2 + AE'E)q'q if - 2q' 2 q 2 if 

+ [4q' 2 q 2 + (q' 2 + q 2 )(E'E - M 2 )} if + 2(E'E + M 2 )q'q if} 
55 V t { = C tt M{[E'(q' 2 + q 2 ) + E(3q' 2 - q 2 )] if - 2(E' + E)q'q if} 
™V t t = C tt M{[E(q' 2 + q 2 ) + E'(3q 2 - q' 2 )] if - 2(E' + E)q'q if} 



with 



v 

Vector-tensor coupling 




Gtt -^Ml ™\A-iF- ( B41 ) 





— C v t 


M[(q' 2 + q 2 ) if ~ 


2q'q if] 


Xi 


— C vt 


M[-{q' 2 + q 2 ) if 


+ Zq'q if 


12 Ki 


— C vt 


M[Qq'q if - 3(q' 2 


+ I 2 ) If] 


3 Xi 


— C v t 


M[2q'q if - (q' 2 ^ 


V)if] 


55j/J 
v vt 


— C v t 


(E'q 2 + 3Eq' 2 ) if 




66t/J 


= Cut 


(Eq' 2 + 3E'q 2 ) if 






(B42) 



with 

Cvt = ^£i£-U—LlJ — J—. (B43) 
4vrM p 2vrM 2 \ E'\ E K ' 

Note that in the p potential, M p is a scaling mass associated with the tensor-coupling 
constant f p . For this scaling mass, the same is to be used in pp, np, and nn scattering. 

More common in nuclear physics is the representation of two-nucleon states in terms of an 
\LSJM) basis, where S denotes the total spin, L the total orbital angular momentum, and 
J the total angular momentum with projection M. In this basis, we denote the potential 
by Vl> S l = {L' SJM\V\LSJM) . These are obtained from the above helicity state matrix 
elements by the following unitary transformation: 
Spin singlet 

yJO = OyJ _ ( B44 ) 
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Uncoupled spin triplet 



V 



ji 
j, j 



lyj 



(B45) 



Coupled triplet states 
1 



V 



ji 



yJl 

V J+1,.J+1 



yJl 

V J-1,J+1 



yJl 

V J+1,J-1 



2J + 1 
1 

2J + 1 
1 

2J + 1 
1 



j 12 v J + (J + 1) 3 V J + p{j + 1)( 5 V J + m v J ) 
(j + 1) 1 2 v J + j u v J - v /j(j + i)( 5 V J + 6f V J } 



j( j + 1)( X V J - u v J ) - j 55 v J + ( j + l) 6 V J ) 
J(j + i)( X V J - 3 V J ) + (J + l) 5 V J - J m v J ) 



2J+1 

The final charge-dependent potentials are 

7(JViiV 2 ) = ^^(A^) + X! ^(A^)] 



(B46) 



(B47) 



with N1N2 either pp, nn, or np. The nucleon mass referred to by M(N\N2) in the above 
equation is fixed as follows 



Af (pp) = M p 
Af (nn) = M n 
M(np) = M e 



M p M n 



938.91875 MeV 



(B48) 
(B49) 

(B50) 



with the precise values for M p and M n given in Table |. The charge-dependent OPE poten- 
tials are given by 



V OPE (pp) = VMM p ),m n0 ,M p ] 
V OPE (nn) = ^(MJ^.MJ 



V OPE (np,T = 1) = —V„ [g-K (Af ) , m^o , M] + 214 [#tt (Af ) , m v ± , Af ] 



(B51) 
(B52) 
(B53) 
(B54) 



V OPE (n Pl T = 0) = —Vn [g w ( Af ) , m n o , Af] - 2 V. [^(Af), m ff±J Af] , 
with m T o and m n ± as given in Table |. Most modern determinations |T^| of the nNN 



coupling constant yield a value for the so-called pseudovector coupling constant, f n [28 



Assuming that f n is fundamentally constant, then g n has a small charge dependence, since 
the two coupling constants are related by 



gl(M) 4M 2 f 2 



An 



m 2 ± An 



(B55) 



with M the mean of the masses of the two nucleons involved in the nNN vertex. We take 
this very small effect into account by using in our V OPE the nNN coupling constant 



gl(M) _ M 2 gl 



An 



Ml An 



(B56) 
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with 



47T 

Defining, 

q 2 4M„ 2 f 2 
4/T m 2 ± 47r 



■ h -13.6. (B57) 



(B58) 



recovers Eq. (|B55| ). 

Since we use units such that h — c — 1, energies, masses and momenta are in units 
of MeV. The potential is in units of MeV" 2 . The conversion factor is he = 197.327053 
MeV fm. If the user wants to relate our units and conventions to the ones used by other 
practitioners, he/she should compare our Eq. ( |A25| ) and our phase shift relation Eq. ( |A33 ) 



with the corresponding equations used by others. A FORTRAN77 computer code for the 
CD-Bonn potential is available from the author. 



APPENDIX C: POTENTIAL PARAMETERS 

For the 'basic' mesons, ir, u, and p, we use, in general, the parameters shown in Table |. 
Note that (except for the cutoff masses) these parameters are determined from empirical 
or semi-empirical sources and, therefore, they are not free parameters of our model. The 
intermediate range attraction is described by two scalar isoscalar bosons, o\ and 02, that 
are also used for the fine-tuning of individual partial waves. The a parameters are given in 
Table [XVI| for the pp T = 1 potential and in Table |XVII| for the T = np potential. For 



partial waves with J > 6, we take gljAn = 2.3 and m ai = 452 MeV. The cutoff mass for 
the two a is A CTl = A CT2 = 2.5 GeV, in all partial waves. In two cases, we vary the cutoff 
parameter of one of the 'basic' mesons: in X P\ we apply A w — > 00 (i. e., the u cutoff is 
omitted), and in 3 P2/ 3 F 2 we use A„. = 3.0 GeV; otherwise, the same cutoff masses (namely, 
the ones shown in Table | and A ai = A CT2 = 2.5 GeV) are used in all partial waves. 

The nn T = 1 potential is constructed by starting from the pp T = 1 potential, replacing 
the proton mass by the neutron mass and adjusting the coupling constants of the two a 
such that the CSB phase shift differences listed in the last column ('Total') of Table [III 
are reproduced. Thus, the a coupling constants of the nn potential (which are given in 
Table pCVlllp are not free parameters. The procedure for the T = 1 np potential is similar. 
We start from the pp T = 1 potential, replace the proton mass by the average mass given 
in Eq. ( B50|) , apply the appropriate OPE potential [i. e., we replace Eq. ( B51|) by ([B53|) 1, 



and then adjust the a coupling constants such that the CIB phase shift differences listed 
in column 'Total' of Table |V| are reproduced which, again, does not generate any free 
parameters. The exception is the 1 So state where the a parameters are used to minimized 
the x 2 m regard to the np data. The charge- dependence caused by the Bonn Full Model 
produces also a small charge-dependent tensor force of 2it range that can be simulated with 
the help of the p coupling. A noticeable effect occurs only in the coupled 3 P2/ 3 -^2 states where 
we use g 2 /4:7T = 0.844 for nn and g 2 /4:7c = 0.862 for np (in all other states g 2 /4:7c = 0.84). 
Again, these choices are made to reproduce the CSB and CIB predicted by other sources 
and, thus, does not introduce new parameters. 
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The free ('fit') parameters of our model are the ones given in Table [XVI| and [XVI ! 
plus two parameters for 1 S np and the cutoff masses which adds up to a total of 43 free 
parameters. 



APPENDIX D: DEUTERON CALCULATIONS 

In momentum space, the deuteron wave function is given by 

*f (k) = [V>o(A0Xi M (k) + UWli^)] C (Dl) 

where yj^ik) are the normalized eigenfunctions of the two-nucleon orbital angular mo- 
mentum L, spin S, and total angular momentum J with projection M; C^ T denotes the 
normalized eigenstates of the total isospin T with projection Mt of the two nucleons. The 
normalization is 



dkk 2 



1. 



(D2) 



The wave functions are obtained by solving the bound state equation which is the ho- 
mogenous version of the scattering equation ( |A18| ): 



V>(k) 



M 



-r 



k 2 



d 3 k'V{k, k')^(k') . 



(D3) 



Note that the deuteron is a pole in the 5* matrix at q — ij. Since we use relativistic 
kinematics in np scattering [cf. Eq. ( A40 )], consistency requires that we determine 7 based 
upon relativistic kinematics which is, 



M d = M p + M n -B d = ^M 2 - 7 2 + \jMl - 7 2 



(D4) 



where M d denotes the deuteron rest mass and Bd the binding energy. The formal solution 
of Eq. (0) is 



7 



AM'M 2 - (M 2 d - M 2 p - M z n Y /AM, 



(D5) 



and, using Bd = 2.224575 MeV and he = 197.327053 MeV fm, the accurate numerical value 
for 7 comes out to be 



7 = 0.2315380 fm" 1 . 
To obtain more insight into 7 2 , we rewrite Eq. (P5|) in factorized form, 



(D6) 



4M d V= (M n + M 1 



M' 



Ml - (M n - M p f 



= B d (4M-B d )(M 2 d -6M 2 ) 
where we introduce the average nucleon mass, 



— _ M p + M n 



938.91897 MeV, 



(D7) 



(D8) 
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and the nucleoli mass difference 5M = M n - M p = 1.29332 MeV, and used M d = 2M - B d 
From this we get 



7 



MB d 1 - 



B d \ 



AM J 



1 



5M 2 



(D9) 



and, in terms of twice the reduced nucleon mass, M, which is defined by 

SM 2 ^ 



- 2 .\ /.,.!/„ 
M= 7 



M„ + M n 



we finally obtain 



7 



M l- 

MB d ( 1 - 
MB d 



938.91852 MeV 



(DIO) 



B d \ 


1 - 


8M 2 
~Mj 


4MJ 


1 - 


SM 2 








B d \ 






4MJ 


. 





(Dll) 

The approximation involved in Eq. ( pil[ ) is good to one part in 10 9 . Therefore, this equation 
reproduces the exact value for 7 to all digits given in Eq. ( |D6|) . One can now identify the 
term MB d as the non-relativistic approximation to 7 2 and the factor (1 — B d /4M) as the 
essential relativistic correction. 

Partial wave decomposition of Eq. (|D3| ) yields for the coupled 3 Si and 3 Di states, 



M 



7 2 + k 2 Jo 
M r 
7 2 + k 2 Jo 



dk'k' 2 [Voo(k, k')Mk') + V 02 (k, k')Mk')} 



dk'k' 2 [V 20 (k, k')Mk') + V 22 (k, k')Mk')] , 



(D12) 



from which ipo and ip 2 are obtained. Considering a finite set of discrete arguments for the 
functions on the l.h.s. and using the same set of momenta to discretize the integrals on the 
r.h.s. produces a matrix equation that is solved easily by the matrix-inversion method ||101 



The momentum-space wave funtions can be Fourier transformed into the configuration- 
space wave functions u and w by 



u L (r) 



- / dkk 2 jL (kr)Mk) 



(D13) 



r y 71 Jo 

with Uo{r) = u(r), u 2 (r) = w(r), and the spherical Bessel functions. The normalization 



is 



dr 



u 2 (r) + w 2 (r) 



1. 



(D14) 



The asymptotic behavior of the wave functions for large values of r are 



u\r\ 



w[r 



A D e 



3 3 
+ 



M dry 



(D15) 
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where As and Ad are known as the asymptotic S- and .D-state normalizations, respectively. 
In addition, one defines the "D/S'-state ratio" r\ = A D /A S . Other deuteron parameters of 
interest are the quadrupole moment 



I f 00 r i— 

Qd = ™ / drr 2 w(r) vSu(r) - w(r) 
zU jo L 



the root-mean-square or matter radius 



rd= \{r drr ' 



and the D-state probability 



u 2 (r) + w 2 {r) 



1/2 



(D16) 



(D17) 



P, 



D 



drw 2 (r). 



(D18) 



The predictions by the CD-Bonn potential for the properties of the deuteron are given in 
Table PCV| numerical values for the wave functions are listed in Table [XIX| and plots are 
shown in Figs. |] and |9|. 

In some applications, it is convenient to have the deuteron wave functions in analytic 
form. Therefore, we present here a simple parametrization of the deuteron functions (that 
was first introduced in Ref. [|110|| ). The ansatz for the analytic version of the r-space wave 
functions is 



u «( r ) = E C i ex P(~ m i r ) 

3=1 
n 

w a{r) = D 3 exp(-m j r) 
3=1 



3 3 
1 + + 



rrijr (rajr)' 



The corresponding momentum space wave functions are 



(D19) 
(D20) 



~{q z + raj 



,1/2^ ^3 



j=i y " l j 



(D21) 
(D22) 



The boundary conditions u a (r) —>■ r and w a (r) — > r 3 as r — > lead to one constraint for the 
Cj and three constraints for the Dj ||110|| , namely 



C n 



n-1 

EQ 



(D23) 
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2 \ ^ ^3 
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(D24) 
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and two other relations obtained by circular permutation of n — 2, n — 1, n. The masses are 



mj = 7 + (J - l)mo (D25) 

with m = 0.9 fm -1 and 7 given in Eq. ( |136| ) . The parameters are given in Table |XX| . The 
constraints Eqs. (P23|) and (P24|) must be enforced by double precision (i. e., to about 15 
decimal digits), otherwise the wave function is not reproduced correctly for r < 0.5 fm. This 
applies, particularly, to the D wave. The accuracy of the parametrization is characterized 



by 

|y dr [u(r) -u a (r)f\ = 2.2 x 10" 4 (D26) 



and 



|^°° dr [w(r) - w a (r)f = 1.1 x 10" 4 . (D27) 

Data files for the deuteron wave functions in r-space as well as in momentum space can be 
obtained from the author upon request. 
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TABLES 



TABLE I. Basic constants and parameters adopted for the CD-Bonn potential. 

Particle Mass (MeV) g 2 /4^ f/g A (GeV) 

tt ± 139.56995 13.6 1 1.72 

vr° 134.9764 13.6 1.72 

p ± ,p° 769.9 0.84 6.1 1.31 

u 781.94 20.0 0.0 1.5 

Proton (p) 938.27231 

Neutron (n) 939.56563 



TABLE II. Differences between the pp and nn 1 Sq effective range parameters [as defined in 
Eqs. ( |3.4j ) and (3.5)] due to the impact of nucleon mass splitting on the kinetic energy (Kin. en.), 
one-boson exchange diagrams (OBE), and two-boson exchanges (TBE). Total denotes the sum of 
the three contributions and empirical information is given in the last column. 





Kin. en. 


OBE 


TBE 


Total 


Empirical 


AacsB (fin) 


0.263 


-0.030 


1.275 


1.508 


1.6 ± 0.6 


Arcss (fm) 


0.004 


0.000 


0.022 


0.026 


0.10 ± 0.12 
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TABLE III. 


Difference 5 nn — 5 PP (in 


degrees) due to the impact of nucleon mass splitting on 


kinematics [30 , 


one-boson exchange dia£ 


5 i ams tvJ-D-Cji, alio iwo-uoson 


excnanges 


[ LE>lh). J-Oiai IS 


the sum of all. 










Ti ab (MeV) 


Kinematics 


OBE 


TBE 


Total 






l S 






0.38254 


0.404 


-0.045 


1.795 


2.154 


1 


0.324 


-0.036 


1.440 


1.728 


5 


0.165 


-0.018 


0.785 


0.932 


10 


0.114 


-0.013 


0.591 


0.692 


25 


0.062 


-0.006 


0.408 


0.464 


50 


0.031 


-0.001 


0.310 


0.340 


100 


0.003 


0.005 


0.239 


0.247 


150 


-0.013 


0.010 


0.206 


0.203 


200 


-0.023 


0.014 


0.185 


0.176 


300 


-0.039 


0.021 


0.160 


0.142 












5 


0.006 


0.001 


0.001 


0.008 


10 


0.013 


0.003 


0.002 


0.018 


25 


0.022 


0.010 


0.008 


0.040 


50 


0.021 


0.021 


0.014 


0.056 


100 


0.004 


0.036 


0.020 


0.060 


150 


-0.011 


0.045 


0.024 


0.058 


200 


-0.022 


0.052 


0.024 


0.054 


300 


-0.040 


0.063 


0.025 


0.048 












5 


-0.003 


0.000 


0.002 


-0.001 


10 


-0.006 


0.000 


0.004 


-0.002 


25 


-0.011 


0.001 


0.012 


0.002 


50 


-0.017 


0.002 


0.027 


0.012 


100 


-0.026 


0.006 


0.049 


0.029 


150 


-0.033 


0.009 


0.065 


0.041 


200 


-0.039 


0.011 


0.076 


0.048 


300 


-0.050 


0.016 


0.090 


0.056 












10 


0.001 


0.000 


0.000 


0.001 


25 


0.002 


0.000 


0.002 


0.004 


50 


0.005 


0.000 


0.006 


0.011 


100 


0.011 


0.002 


0.019 


0.032 


150 


0.016 


0.005 


0.033 


0.054 


200 


0.019 


0.010 


0.046 


0.075 
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300 


0.022 


0.022 


0.068 


0.112 












5 


0.001 


0.000 


0.001 


0.002 


10 


0.003 


0.000 


0.004 


0.007 


25 


0.010 


0.001 


0.013 


0.024 


50 


0.021 


0.002 


0.031 


0.054 


100 


0.032 


0.006 


0.062 


0.100 


150 


0.036 


0.010 


0.081 


0.127 


200 


0.035 


0.015 


0.093 


0.143 


300 


0.032 


0.023 


0.105 


0.160 
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TABLE IV. Differences between the pp and np 1 Sq effective range parameters [as denned in 



Eqs. ( |3.6| ) and (3.7)] produced by various CIB mechanisms and phenomenology. Total is the sum 
of all contributions listed left of column 'Total'. AM denotes all effects caused by nucleon mass 
splitting. Empirical information is given in the last column. 





AM 


OPE 


TBE 


7T7 


phenom. 


Total 


Empirical 


Aa C iB (fm) 


0.754 


3.035 


1.339 


-0.405 


1.555 


6.278 


6.44 ± 0.40 


Arc/s (fin) 


0.013 


0.092 


0.016 


-0.004 


0.057 


0.174 


0.08 ± 0.06 



TABLE V. Difference 

$np ~ 3pp (in degrees) in the Sq state as produced by various CIB 



mechanisms and phenomenology. AM stands for all effects caused by nucleon mass splitting [41] 



Total is the sum of all contributions listed left of column 'Total'. 'All' denotes the sum of Total 



and Coulomb, 


where Coulomb is 


the difference 5 PP — 


5 C . 

pp 








Ti ab (MeV) 


AM 


OPE 


TBE 




phenom. 


Total 


Coulomb 


All 


0.38254 


1.077 


3.541 


1.655 


-0.412 


1.953 


7.814 


32.085 


39.894 


1 


0.859 


2.851 


1.260 


-0.305 


1.521 


6.186 


23.114 


29.300 


5 


0.468 


1.650 


0.654 


-0.152 


0.982 


3.602 


5.219 


8.821 


10 


0.350 


1.271 


0.482 


-0.106 


0.909 


2.906 


1.896 


4.802 


25 


0.240 


0.875 


0.320 


-0.058 


0.970 


2.347 


-0.044 


2.304 


50 


0.182 


0.656 


0.233 


-0.028 


1.142 


2.185 


-0.589 


1.597 


100 


0.139 


0.513 


0.165 


-0.002 


1.433 


2.248 


-0.772 


1.476 


150 


0.119 


0.469 


0.130 


0.012 


1.656 


2.386 


-0.796 


1.590 


200 


0.108 


0.457 


0.103 


0.021 


1.839 


2.528 


-0.796 


1.733 


300 


0.094 


0.477 


0.058 


0.034 


2.124 


2.787 


-0.782 


2.005 
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TABLE VI. Difference $np — $pp (in degrees) for partial waves with L > as produced by 

various CIB mechanisms. Notation as in Table | V|. 

Ti ab (MeV) AM OPE TBE vr 7 Total Coulomb All 



1 


0.000 


-0.030 


0.000 


0.000 


-0.030 


0.073 


0.043 


5 


0.000 


-0.230 


-0.003 


0.000 


-0.233 


0.262 


0.029 


10 


0.000 


-0.448 


-0.009 


0.000 


-0.457 


0.353 


-0.104 


25 


0.012 


-0.770 


-0.027 


-0.017 


-0.802 


0.320 


-0.481 


50 


0.032 


-0.846 


-0.050 


-0.050 


-0.914 


0.111 


-0.803 


100 


0.050 


-0.742 


-0.074 


-0.087 


-0.853 


-0.142 


-0.996 


150 


0.050 


-0.649 


-0.083 


-0.104 


-0.786 


-0.255 


-1.041 


200 


0.047 


-0.586 


-0.088 


-0.113 


-0.740 


-0.314 


-1.054 


300 


0.045 


-0.513 


-0.096 


-0.125 


-0.689 


-0.369 


-1.058 


1 


0.000 


0.016 


"Pi 

0.000 


0.000 


0.016 


-0.043 


-0.026 


5 


0.002 


0.110 


0.001 


-0.002 


0.111 


-0.140 


-0.028 


10 


0.004 


0.193 


0.003 


-0.002 


0.198 


-0.187 


0.011 


25 


0.006 


0.298 


0.008 


0.003 


0.315 


-0.224 


0.091 


50 


0.008 


0.330 


0.018 


0.016 


0.372 


-0.240 


0.133 


100 


0.016 


0.307 


0.038 


0.038 


0.399 


-0.265 


0.133 


150 


0.022 


0.274 


0.055 


0.054 


0.405 


-0.287 


0.118 


200 


0.028 


0.246 


0.069 


0.064 


0.407 


-0.303 


0.103 


300 


0.033 


0.202 


0.099 


0.077 


0.411 


-0.325 


0.085 


5 


0.000 


-0.009 


0.000 


0.000 


-0.009 


0.007 


-0.002 


10 


0.000 


-0.024 


0.000 


0.000 


-0.024 


0.015 


-0.009 


25 


0.000 


-0.049 


0.001 


0.001 


-0.047 


0.031 


-0.016 


50 


0.002 


-0.043 


0.005 


-0.002 


-0.038 


0.049 


0.011 


100 


0.014 


0.003 


0.013 


-0.011 


0.019 


0.071 


0.090 


150 


0.024 


0.041 


0.023 


-0.018 


0.070 


0.081 


0.151 


200 


0.034 


0.068 


0.030 


-0.025 


0.107 


0.083 


0.190 


300 


0.045 


0.095 


0.042 


-0.033 


0.149 


0.073 


0.222 


5 


0.000 


-0.009 


-0.001 


0.000 


-0.010 


0.049 


0.040 


10 


0.001 


-0.028 


-0.002 


0.000 


-0.029 


0.094 


0.065 


25 


0.004 


-0.090 


-0.005 


-0.001 


-0.092 


0.188 


0.097 


50 


0.017 


-0.162 


-0.011 


-0.006 


-0.162 


0.257 


0.095 


100 


0.043 


-0.211 


-0.024 


-0.020 


-0.212 


0.260 


0.048 


150 


0.058 


-0.210 


-0.032 


-0.030 


-0.214 


0.221 


0.007 


200 


0.065 


-0.196 


-0.035 


-0.037 


-0.203 


0.184 


-0.019 
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300 0.072 -0.169 -0.034 -0.044 -0.175 0.130 -0.044 



3 F 2 

10 0.000 -0.004 0.000 0.000 -0.004 0.001 -0.002 

25 0.000 -0.019 0.000 0.000 -0.019 0.004 -0.015 

50 0.000 -0.043 0.000 0.001 -0.042 0.007 -0.036 

100 0.000 -0.068 0.000 0.002 -0.066 0.008 -0.058 

150 0.003 -0.081 -0.001 0.002 -0.077 0.007 -0.070 

200 0.007 -0.090 -0.001 0.002 -0.082 0.003 -0.079 

300 0.008 -0.099 -0.001 0.002 -0.090 -0.009 -0.098 



(2 



5 


0.000 


0.011 


0.000 


0.000 


0.011 


-0.008 


0.004 


10 


0.001 


0.034 


0.000 


-0.001 


0.034 


-0.016 


0.018 


25 


0.002 


0.086 


0.000 


-0.002 


0.086 


-0.028 


0.058 


50 


-0.001 


0.111 


0.003 


0.001 


0.114 


-0.025 


0.089 


100 


-0.004 


0.087 


0.007 


0.010 


0.100 


-0.003 


0.097 


150 


-0.004 


0.051 


0.010 


0.018 


0.075 


0.017 


0.092 


200 


-0.001 


0.020 


0.012 


0.024 


0.055 


0.032 


0.087 


300 


0.008 


-0.020 


0.014 


0.032 


0.034 


0.047 


0.080 
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TABLE VII. pp phase shifts in degrees. 



T lab (MeV) 


'So 


3 Po 


3 Pi 


l D 2 


3 P 2 


3 F2 


£2 


3 ^3 


'G 4 


3 F 4 


1 


32.79 


0.13 


-0.08 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


5 


54.85 


1.58 


-0.90 


0.04 


0.22 


0.00 


-0.05 


0.00 


0.00 


0.00 


10 


55.20 


3.72 


-2.05 


0.17 


0.66 


0.01 


-0.20 


-0.03 


0.00 


0.00 


25 


48.63 


8.58 


-4.90 


0.70 


2.50 


0.10 


-0.81 


-0.23 


0.04 


0.02 


50 


38.86 


11.54 


-8.31 


1.71 


5.84 


0.33 


-1.73 


-0.70 


0.15 


0.12 


100 


24.91 


9.57 


-13.37 


3.77 


10.97 


0.78 


-2.72 


-1.53 


0.42 


0.50 


150 


14.73 


4.76 


-17.62 


5.67 


13.98 


1.10 


-2.99 


-2.12 


0.69 


1.04 


200 


6.58 


-0.49 


-21.49 


7.26 


15.68 


1.27 


-2.88 


-2.48 


0.97 


1.63 


250 


-0.29 


-5.62 


-25.05 


8.55 


16.63 


1.26 


-2.59 


-2.68 


1.26 


2.19 


300 


-6.26 


-10.48 


-28.36 


9.54 


17.12 


1.08 


-2.21 


-2.75 


1.55 


2.69 


350 


-11.56 


-15.04 


-31.45 


10.27 


17.33 


0.73 


-1.80 


-2.72 


1.83 


3.11 



TABLE VIII. nn phase shifts in degrees. 



T iab (MeV) 


'So 


3 p 

Pq 


3 Pi 


l D 2 


3 P2 


3 F2 


£2 


3 ^3 


'G 4 


3 Fa 


1 


57.63 


0.21 


-0.12 


0.00 


0.02 


0.00 


0.00 


0.00 


0.00 


0.00 


5 


61.00 


1.85 


-1.04 


0.05 


0.27 


0.00 


-0.06 


-0.01 


0.00 


0.00 


10 


57.79 


4.10 


-2.24 


0.18 


0.76 


0.01 


-0.22 


-0.04 


0.00 


0.00 


25 


49.05 


8.94 


-5.13 


0.74 


2.71 


0.11 


-0.85 


-0.24 


0.04 


0.02 


50 


38.61 


11.71 


-8.54 


1.77 


6.15 


0.34 


-1.76 


-0.71 


0.16 


0.12 


100 


24.38 


9.49 


-13.60 


3.88 


11.33 


0.79 


-2.73 


-1.55 


0.42 


0.52 


150 


14.14 


4.56 


-17.87 


5.80 


14.32 


1.11 


-2.97 


-2.13 


0.70 


1.06 


200 


5.96 


-0.75 


-21.74 


7.42 


16.01 


1.28 


-2.85 


-2.49 


0.98 


1.66 


250 


-0.92 


-5.92 


-25.31 


8.72 


16.94 


1.27 


-2.54 


-2.68 


1.28 


2.23 


300 


-6.90 


-10.80 


-28.63 


9.72 


17.42 


1.09 


-2.15 


-2.74 


1.57 


2.73 


350 


-12.21 


-15.38 


-31.72 


10.46 


17.60 


0.73 


-1.74 


-2.70 


1.86 


3.15 
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TABLE IX. T = 1 np phase shifts in degrees. 



T lab (MeV) 


'So 


3 ^o 


3 Pi 


l D 2 


3 P 2 


3 F 2 


£2 


3 ^3 


l G 4 


3 F 4 


1 


62.09 


0.18 


-0.11 


0.00 


0.02 


0.00 


0.00 


0.00 


0.00 


0.00 


5 


63.67 


1.61 


-0.93 


0.04 


0.26 


0.00 


-0.05 


0.00 


0.00 


0.00 


10 


60.01 


3.62 


-2.04 


0.16 


0.72 


0.01 


-0.18 


-0.03 


0.00 


0.00 


25 


50.93 


8.10 


-4.81 


0.69 


2.60 


0.09 


-0.76 


-0.20 


0.03 


0.02 


50 


40.45 


10.74 


-8.18 


1.73 


5.93 


0.30 


-1.64 


-0.62 


0.13 


0.11 


100 


26.38 


8.57 


-13.23 


3.86 


11.01 


0.72 


-2.63 


-1.42 


0.39 


0.48 


150 


16.32 


3.72 


-17.51 


5.82 


13.98 


1.03 


-2.90 


-1.98 


0.67 


1.01 


200 


8.31 


-1.55 


-21.38 


7.45 


15.66 


1.19 


-2.79 


-2.33 


0.96 


1.59 


250 


1.59 


-6.68 


-24.96 


8.76 


16.59 


1.17 


-2.50 


-2.51 


1.26 


2.15 


300 


-4.25 


-11.54 


-28.27 


9.76 


17.08 


0.98 


-2.13 


-2.57 


1.56 


2.65 


350 


-9.44 


-16.10 


-31.37 


10.49 


17.28 


0.62 


-1.72 


-2.53 


1.85 


3.06 



TABLE X. T = np phase shifts in degrees. 



T lab (MeV) 


'Pi 


3 s 1 


3 Di 


ei 


3 D 2 


^3 


3 ^3 


3 G 3 


£3 


3 G 4 


1 


-0.19 


147.75 


-0.01 


0.11 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


5 


-1.49 


118.18 


-0.18 


0.68 


0.22 


-0.01 


0.00 


0.00 


0.01 


0.00 


10 


-3.05 


102.62 


-0.68 


1.17 


0.85 


-0.07 


0.01 


0.00 


0.08 


0.01 


25 


-6.35 


80.63 


-2.80 


1.81 


3.72 


-0.42 


0.05 


-0.05 


0.55 


0.17 


50 


-9.73 


62.73 


-6.44 


2.13 


8.97 


-1.11 


0.33 


-0.26 


1.61 


0.72 


100 


-14.43 


43.06 


-12.25 


2.45 


17.22 


-2.15 


1.45 


-0.94 


3.49 


2.17 


150 


-18.33 


30.47 


-16.50 


2.79 


22.09 


-2.87 


2.70 


-1.76 


4.83 


3.64 


200 


-21.77 


20.95 


-19.68 


3.18 


24.51 


-3.48 


3.70 


-2.60 


5.76 


4.99 


250 


-24.84 


13.21 


-22.12 


3.60 


25.36 


-4.08 


4.31 


-3.39 


6.40 


6.18 


300 


-27.57 


6.65 


-24.03 


4.00 


25.21 


-4.73 


4.54 


-4.09 


6.83 


7.21 


350 


-30.00 


0.92 


-25.53 


4.38 


24.44 


-5.45 


4.44 


-4.71 


7.14 


8.07 
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TABLE XI. Scattering lengths (a) and effective ranges (r) in units of fm. 



CD-Bonn 



Experiment 



Reference(s) 







-7.8154 


-7.8149 ± 0.0029 


2.773 


2.769 ± 0.014 


-17.4602 




2.845 




-18.9680 


-18.9 ± 0.4 


2.819 


2.75 ± 0.11 


-23.7380 


-23.740 ± 0.020 


2.671 


( 2.77 ± 0.05 ) 



a c 
r c 
pp 
a N 

r N 
PP 
N 
nn 
N 
nn 



up 
r np 



m 

[27| 



(H 

n 



5.4196 
1.751 



3 S! 

5.419 ± 0.007 
1.753 ± 0.008 
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TABLE XII. After-1992 pp data below 350 MeV included in the 1999 pp database. Error refers 
to the normalization error. This table contains 1113 observables and 32 normalizations resulting 
in a total of 1145 data. 



Ti ab (MeV) 



# observable 



Error (%) 



Institution(s) 



Reference 



25.68 
25.68 
25.68 
197.4 
197.4 
197.4 
197.4 
197.4 
197.8 
197.8 
197.8 
197.8 
197.8 
197.8 
197.8 
197.8 
197.8 
250.0 
250.0 
250.0 
250.0 
280.0 
280.0 
280.0 
280.0 
294.4 
294.4 
294.4 
294.4 
310.0 
310.0 
310.0 
310.0 
350.0 
350.0 
350.0 
350.0 



8 D 
6 R 
2 A 
41 P 
41 A xx 

41 Ayy 

41 A xz 

39 A zz 
14 P 

14 A xx 

14 Ayy 

14 A xz 
10 D 
5 R 
5 R' 
5 A 
5 A' 
41 P 

41 A xx 

41 Ayy 

41 A xz 
41 P 
41 A xx 

41 Ayy 

41 

40 P 

40 i ra 

40 Ayy 

40 Am 
40 P 
40 A xx 
40 A w 
40 A xz 
40 P 
40 

40 Ayy 

40 ^ 



1.3 
1.3 
1.3 
1.3 
2.5 
2.5 
2.5 
2.0 
1.3 
2.4 
2.4 
2.4 
None 
None 
None 
None 
None 
1.3 
2.5 
2.5 
2.5 
1.3 
2.5 
2.5 
2.5 
1.3 
2.5 
2.5 
2.5 
1.3 
2.5 
2.5 
2.5 
1.3 
2.5 
2.5 
2.5 



Erlangen, Zurich, PSI 
Erlangen, Zurich, PSI 
Erlangen, Zurich, PSI 
onsin, 1 
□nsin, 1 
□nsin, 1 
onsin, 1 
onsin, 1 
onsin, 1 
onsin, 1 
onsin, 1 
onsin, 1 
IUCF 
IUCF 
IUCF 
IUCF 
IUCF 



Wisconsin, 


IUCF 


Wisconsin, 


IUCF 


Wisconsin, 


IUCF 


Wisconsin, 


IUCF 


Wisconsin, 


IUCF 


Wisconsin, 


IUCF 


Wisconsin, 


IUCF 


Wisconsin, 


IUCF 


Wisconsin, 


IUCF 



IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 


IUCF, 


Wisconsin 



m 
m 
m 



m 
m 

B 

m 



m 
m 
\m 
m 
m 

\B 

m 

m 
m 

\m 



m 

m 
m 
\m 
m 
m 



m 
m 
m 

\m 
m 

\m 
m 
m 
m 
m 
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TABLE XIII. After-1992 np data below 350 MeV included in the 1999 np database. Error refers 
to the normalization error. This table contains 524 observables and 20 normalizations resulting in 
a total of 544 data. 

Ti a b (MeV) # observable Error (%) Institution(s) Reference 



3.65-11.6 9 Aa T None TUNL fH 

4.98-19.7 6 Aa L None TUNL @ 

4.98-17.1 5 Act None TUNL @ 

14.11 6 a 0.7 Tubingen || 

15.8 1 A None Bonn @ 

16.2 1 Ao- r None Prague || 

16.2 1 Aa L None Prague || 

175.26 84 P Float" TRIUMF || 

203.15 100 P 4.7 TRIUMF || 

217.24 100 P 4.5 TRIUMF || 

260.0 8 R t 3.0 PSI @ 

260.0 8 A t 3.0 PSI || 

260.0 3 A t 3.0 PSI f6§ 

260.0 8 A 3.0 PSI H 

260.0 3 A 3.0 PSI || 

260.0 8 P 2.0 PSI |H 

260.0 3 P 2.0 PSI || 

261.00 88 P 4.1 TRIUMF |63j 

312.0 24 P 4.0 SATURNE §| 

312.0 11 A zz 4.0 SATURNE |6j 

318.0 8 R t 3.0 PSI H| 

318.0 8 A t 3.0 PSI || 

318.0 5 A t 3.0 PSI || 

318.0 8 A 3.0 PSI || 

318.0 5 A 3.0 PSI |H 

318.0 8 P 2.0 PSI || 

318.0 5 P 2.0 PSI 



a This data set is floated because all current phase shift analyses and np potentials predict a 
norm that is about 4 standard deviations off the experimental normalization error of 4.9%, 
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TABLE XIV. x 2 /datum for the CD-Bonn potential, the Nijmegen phase shift analysis [42], 
and the Argonne Vis potential |32| in regard to various databases discussed in the text. 





CD-Bonn 


Nijmegen 


Argonne 




potential phase shift analysis 


Vis potential 




proton-proton data 






1992 pp database (1787 data) 


1.00 


1.00 


1.10 


After-1992 pp data (1145 data) 


1.03 


1.24 


1.74 


1999 pp database (2932 data) 


1.01 


1.09 


1.35 




neutron-proton data 






1992 np database (2514 data) 


1.03 


0.99 


1.08 


After-1992 np data (544 data) 


0.99 


0.99 


1.02 


1999 np database (3058 data) 


1.02 


0.99 


1.07 




pp and np data 






1992 NN database (4301 data) 


1.02 


0.99 


1.09 


1999 NN database (5990 data) 


1.02 


1.04 


1.21 
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TABLE XV. Deuteron properties. 



CD-Bonn Empirical 



Reference(s) 



Binding energy (MeV) 
Deuteron effective range = p(—Bd 
Asymptotic S state As (fm -1 / 2 ) 
Asymptotic D/S state 77 
Matter radius (fm) 
Quadrupole moment Qd (fm 2 ) 
D-state probability Pd 



-B d ) (fm) 



2.224575 

1.765 

0.8846 

0.0256 

1.966 

0.270 a 

4.85 



2.224575(9) 

1.765(9) 

0.8846(9) 

0.0256(4) 

1.971(6) 

0.2859(3) 



[23§9p| 



I sap 1 



a Without meson current contributions and relativistic corrections. 
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TABLE XVI. Parameters of the scalar isoscalar bosons, a\ and a 2 , for the pp T = 1 potential. 
An asterix denotes the default which are the 1 Sq parameters. The boson masses m ai and m a2 are 
in units of MeV. 





(/^/4ir (m CTl ) 


9l 2 /4vr (m CT2 ) 


'So 


4.24591 (452) 


17.61 (1225) 




7.866 (560) 


* (*) 


3 P1 


2.303 (424) 


* (*) 


3 P2 


4.166 (470) 


24.80 (*) 


l D 2 


2.225 (400) 


190.7 (*) 


r 2 , -^3 


1.5 (*) 


56.21, 74.44 (793) 


3 F 4 , 3 tf 4 


3.8 (*) 


* (*) 


*G 4 


* (*) 




3 # 5 


* (*) 





TABLE XVII. Parameters of the scalar isoscalar bosons, o\ and a 2 , for the T = np potential. 
An asterix denotes the default which are the 3 Si parameters. The boson masses m ai and m a2 are 



in units of MeV. 




glj4ir (m ffl ) 


d 2 /4vr (m CT2 ) 


3 Si 


0.51673 (350) 


14.01164 (793) 


'Pi, 3 D 2 


0.81, 0.53 (*) 


71.5, 154.5 (1225) 


3 D 1 


0.575 (*) 




3 D 3 


3.4 (452) 




'F 3 


0.73 (*) 




3 G 3 


0.29 (*) 




3 G 4 


0.62 (*) 




3 G 5 , 3 / 5 


0.96 (*) 




'H 5 


* (*) 
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TABLE XVIII. Coupling constants of the scalar isoscalar bosons, o\ and 02, for the T = 1 np 
and nn potentials. Note that these are not free parameters (except for 1 So np). The boson masses 
are the same as for the pp T = 1 potential (Table XVI ) . 



— neutron-proton — — neutron-neutron — 







<4/4vr 


9l t /^ 


<4/4tt 


1 5'o 


3.96451 


22.50007 


4.26338 


17.54 


3 Po 


7.866 


5.8 


7.892 


16.747 


3 Pi 


2.346 


19.22 


2.326 


17.61 


3 P2 


4.194 


24.562 


4.180 


24.737 


l D 2 


2.236 


189.7 


2.241 


190.7 


3 7? 3 771 


1.573, 1.53 


56.21, 74.85 


1.522, 1.53 


56.28, 74.44 




3.8115, 3.85 


17.61 


3.81, 3.83 


17.61 




4.27591 




4.284 




3 H 5 


4.24591 




4.24591 
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TABLE XIX. Deuteron wave functions. 



r (fm) 


u(r) (fm 


1/2) 


w 


(r) (fm-V2) 




r (fm) 




u{r) (fm 


1/2) 


w(r) (fm 


1/2) 


0.1006 


- 01 





.3040616 


-02 


-0, 


.1372766 


- 05 


o. 


.2706 


+ 


01 





,4575506 


+ 00 


0. 


.1072196 + 00 


0.2006; 


- 01 


0. 


.6073136 


- 02 


-0. 


.8952156 


- 05 





.2806 


+ 


01 


0. 


.4488376 


+ 00 


0. 


.1025726? + 00 


0.3006 


- 01 


0. 


.9094446 


- 02 


-0. 


2494956; 


-04 





.2906 


+ 


01 


0. 


.4400646? 


+ 00 


0. 


.9807686 


- 01 


0.4006 


- 01 





.1210486; 


- 01 


-0. 


4923126 


-04 





.3006 


+ 


01 


0. 


.4312756 


+ 00 


0. 


.9374536 


- 01 


0.5006 


- 01 


0. 


.1510656; 


- 01 


-0. 


.8042756 


-04 





.3206 


+ 


01 


0. 


.4137786 


+ 00 


0. 


8559236; 


- 01 


omoE 


- 01 





.1810296 


- 01 


-0. 


1166106 


- 03 





.U0E 


+ 


01 


0. 


.3965526? 


+ 00 


0. 


.7812356 


- 01 


0.7006 


- 01 


0. 


,2109846; 


- 01 


-0. 


1555266? 


- 03 





.360 E 


+ 


01 


0. 


.3797276 


+ 00 


0. 


7131766 


-01 


0.8006 


- 01 


0. 


,2409756 


-01 


-0, 


,1948136; 


- 03 





.380 E 


+ 


01 


0. 


.3633876 


+ 00 


0. 


651366,6 


-01 


0.9006 


- 01 


0. 


.2710506; 


- 01 


-0. 


.2320586; 


- 03 


o. 


A00E 


+ 


01 


0. 


.3475836? 


+ 00 


0. 


5953446 


- 01 


0.1006 


+ 00 


0. 


.3012556; 


- 01 


-0. 


.2648716 


- 03 





A20E 


+ 


01 


0. 


.3323436? 


+ 00 


0. 


5446236; 


- 01 


0.2006; 


+ 00 





.6201936 


-01 


0. 


.1556436; 


-03 





.4406 


+ 


01 


0. 


.3176786 


+ 00 


0. 


4987216 


- 01 


0.3006 


+ 00 





.9938766 


-01 


0. 


3350716 


-02 





.4606 


+ 


01 


0. 


.3035926? 


+ 00 


0. 


4571786 


-01 


0.4006 


+ 00 


0. 


.143869£ 


+ 00 


0. 


.1089366 


- 01 


0. 


.4806 


+ 


01 


0. 


.2900786 


+ 00 


0. 


4195656; 


-01 


0.5006 


+ 00 


0. 


.1945456 


+ 00 


0. 


.2355746 


- 01 





.5006 


+ 


01 


0. 


.2771266? 


+ 00 


0. 


.3854876 


-01 


omoE 


+ 00 





.2484546 


+ 00 


0. 


4090686? 


- 01 





.5206 


+ 


01 





.2647216 


+ 00 


0. 


.3545876 


- 01 


0.7006 


+ 00 





.3018416 


+ 00 


0. 


6128086? 


- 01 





.5406 


+ 


01 





.2528496? 


+ 00 


0. 


.3265406; 


- 01 


0.8006 


+ 00 





.3513746 


+ 00 


0. 


8240336? 


- 01 


0. 


.5606 


+ 


01 





.2414916 


+ 00 


0. 


3010566; 


- 01 


0.9006 


+ 00 





.394806,6 


+ 00 


0, 


.1021766 


+ 00 


0. 


.5806 


+ 


01 





,2306296? 


+ 00 


0. 


.2778746 


-01 


0.1006? 


+ 01 





,4310726 


+ 00 


0. 


1191656; 


+ 00 


0. 


moE 


+ 


01 





.2202456 


+ 00 


0. 


2567616 


-01 


O.UOE 


+ 01 





,460046£ 


+ 00 


0. 


1326836? 


+ 00 


0. 


.6506 


+ 


01 





.1962526? 


+ 00 


0. 


.2117176 


-01 


0.1206; 


+ 01 





.4822136; 


+ 00 


0. 


.1426336? 


+ 00 


0. 


.100E 


+ 


01 





.1748466? 


+ 00 


0. 


1756766 


- 01 


0.1306; 


+ 01 





,4983706 


+ 00 


0. 


.1492856? 


+ 00 


0. 


.7506? 


+ 


01 





.1557596? 


+ 00 


0. 


.1466166 


- 01 


0.1406; 


+ 01 





,5094156 


+ 00 


0. 


.153089,6 


+ 00 


0. 


moE 


+ 


01 





,138747.6 


+ 00 


0. 


1230106 


-01 


0.1506; 


+ 01 





.5162226; 


+ 00 


0, 


.1545456 


+ 00 


0. 


.8506 


+ 


01 





,1235896? 


+ 00 


0. 


103699,6 


-01 


0.1606; 


+ 01 


0. 


.5195796 


+ 00 


0. 


1541366? 


+ 00 





.9006 


+ 


01 


0. 


.1100846? 


+ 00 


0. 


.8779936; 


- 02 


0.1706; 


+ 01 


0. 


.5201586; 


+ 00 


0. 


.1522876 


+ 00 





.9506 


+ 


01 


0. 


.9805256 


- 01 


0. 


.7462816 


- 02 


0.1806; 


+ 01 


0. 


.5185246; 


+ 00 


0. 


1493566? 


+ 00 





.1006 


+ 


02 


0. 


.8733546 


- 01 


0. 


.6365656; 


- 02 


0.1906; 


+ 01 


0. 


.5151386 


+ 00 


0. 


.1456386; 


+ 00 





.1056? 


+ 


02 


0. 


.7778916 


-01 


0. 


.5447056; 


- 02 


0.2006 


+ 01 





.5103746 


+ 00 


0. 


,1413676; 


+ 00 





.1106? 


+ 


02 





.6928596; 


-01 


0. 


.4674386; 


- 02 


0.2106 


+ 01 


0. 


.5045336; 


+ 00 


0. 


.1367286? 


+ 00 





.1156 


+ 


02 


0. 


,6171206; 


- 01 


0. 


4021706 


- 02 


0.2206 


+ 01 


0. 


.4978566; 


+ 00 


0, 


.1318646? 


+ 00 





A20E 


+ 


02 


0. 


.5496606? 


- 01 


0. 


3468266? 


- 02 


0.2306 


+ 01 





.4905396; 


+ 00 


0. 


.1268866? 


+ 00 


0. 


.1256? 


+ 


02 


0. 


.4895736 


- 01 


0. 


2997346; 


- 02 


0.2406 


+ 01 


0. 


,482736£ 


+ 00 


0. 


.1218776 


+ 00 





.1306 


+ 


02 


0. 


,4360556; 


- 01 


0. 


2595356 


- 02 


0.2506 


+ 01 


0. 


,4745736 


+ 00 


0. 


.1169016 


+ 00 


0. 


.1356 


+ 


02 


0. 


.388386,6 


- 01 


0. 


2251206; 


-02 


0.2606 


+ 01 


0. 


,4661506 


+ 00 


0, 


.1120046 


+ 00 


0. 


.1406 


+ 


02 


0. 


.3459296 


- 01 


0. 


.1955826 


- 02 
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TABLE XX. Coefficients for the parametrized deuteron wave functions (n = 11). 



3 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 



Cj (fm- 1 ^) 



Dj (fm 



-1/2^ 



0.88472985L> + 00 
-0.26408759L* + 00 
-0.44114404D - 01 
-0.14397512D + 02 

0.85591256D + 02 
-0.31876761D + 03 

0.70336701D + 03 
-0.90049586D + 03 

0.66145441D + 03 
-0.25958894D + 03 
Eq. ([D23D 



0.22623762D - 01 
-0.50471056Z) + 00 

0.562788971? + 00 
-0.16079764D + 02 

0.11126803D + 03 
-0.44667490D + 03 

0.109859071? + 04 
-0.16114995D + 04 



Eq. (D24) 



Eq. (|D 24|) 



Eq. (D24) 
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FIGURES 




FIG. 1. One-boson exchange Feynman diagrams that define the CD-Bonn NN potential. 




i i i i i i i i i 

0.5 1 1.5 2 
q (GeV) 



FIG. 2. Half off-shell 3 £i— 3 Z?i potential. The on-shell momentum q' is held fixed at 265 MeV 
(equivalent to 150 MeV lab. energy), while the off-shell momentum q runs from zero to 2000 MeV. 
The on-shell point (q = 265 MeV) is marked by a solid dot. The solid curve is the relativistic OBE 
amplitude of ir + p exchange. When the relativistic OPE amplitude, Eq. (|2.9|) , is replaced by the 
static/local approximation, Eq. ( |2.10| ), the dashed curve is obtained, and when this approximation 
is also used for the one-p exchange, the dotted curve results. 
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FIG. 3. Differences 5 nn — 5 PP due to the impact of nucleon mass splitting on kinematics (dotted 
line labeled 'Kin'), one-boson exchange diagrams (dashed double-dotted, OBE), and two-boson 
exchanges (dashed, TBE). The solid line ('Tot') represents the total. Notice that each frame has 
a different scale. 
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FIG. 4. Differences $np — 3pp as produced by various CIB mechanisms. Shown are the contri- 
butions from OPE (dashed curve), TBE (dashed-dotted), and irreducible irj exchange (dotted). 
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FIG. 5. The difference 5 np — 5 PP due to the charge-dependence of the strong force (dashed 
curve labeled 'Tot') and (8 PP — 5 pp ) due to the Coulomb force (dotted, Cb). The sum of both is 
represented by the solid line labeled 'All'. 
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FIG. 6. pp phase parameters in partial waves with J < 4. The solid line represents the 
predictions by the CD-Bonn potential. The solid dots and open circles are the results from the 



Nijmegen multi-energy pp phase shift analysis Q| and the VPI single-energy pp analysis SM99 [45], 
respectively. 
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Fig. H continued. 
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FIG. 7. np phase parameters in partial waves with J < 4. The solid line represents the 
predictions by the CD-Bonn potential. The solid dots and open circles are the results from the 
Nijmegen multi-energy np phase shift analysis [42| and the VPI single-energy np analysis SM99 [45], 
respectively. 
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Fig. [7| continued. 
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Fig. [7| continued. 
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FIG. 8. Deuteron wave functions. The family of large curves are u(r) and the family of small 
curves are w(r). The solid lines represent the wave functions generated from the CD-Bonn potential, 
while the dashed and dotted lines are from the Nijmegen-I 31] and Argonne Vis (Hi potentials, 
respectively. 
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FIG. 9. The deuteron wave functions of Fig. || in an alternative representation. The family of 
large curves are u(r)/r and the family of small curves are w(r)/r. 
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